Prof A Kaufmann Ing G Desbazeille METODA DRUMULUI CRITIC Aplicarea metodei PERT şi a variantelor sale la programele de producjie şi de învăjămînt Traducere din limba franceză, după edijia a doua Editura Editura tehnică tehnică Bucureşti — 1971 Traducere şi control ştiinţific: prof Gavrll Sâmboan şi dr loan Tomescu Redactor: Valentina Bucur Tehnoredactor: Vaieriu Morărescu Coperta şi supracoperta: Valentin Vişan La mâthode du chemin critique par A Kaufmann et G Desbazeille © Dunod, Paris, 1969 CUVÎNT ÎNAINTE Iată o lucrare foarte utilă, care vine la timp şi căreia îi dorim o foarte largă popularitate Curiozitatea şi interesul trezite aproape pretutindeni în lume de sistemul PEET sau metoda drumului critic („Criticai Path Scheduling“) în cursul ultimilor ani se bazează pe mai multe cauze Sînt înclinat să cred că succesul ei se datorează în mare parte faptului că ea nu necesită, ca alte metode de cercetare operaţională, o formaţie matematică prealabilă avansată şi că noţiunile la care apelează sînt destul de intuitive De asemenea, ea se asimilează uşor şi, prin urmare, este accesibilă Bineînţeles, ar mai trebui ca această metodă să răspundă unor necesităţi reale şi adesea resimţite de toţi aceia care au de dus la bun sfîrşit realizarea unui proiect sau a unui mare ansamblu Chiar aceasta este efectiv situaţia, deoarece insuficienţele diagramei clasice a lui Gantt au devenit flagrante Lucrarea de faţă are, după părerea mea, un merit esenţial: acela de a fi în primul rînd practică După un prim capitol de recapitulări teoretice asupra teoriei grafurilor şi a problemelor de ordonanţare, autorii abordează noţiuni concrete şi fundamentale, ca antecedenţi şi descendenţi, timpi operatori, interval de fluctuaţie, margini (totale, libere, sigure), pe care practicienii nu vor avea nici o dificultate în a le identifica, referin- du-se la propria lor experienţă Apoi urmează preţioasele exemple numerice datorită cărora algoritmii prezentaţi îşi pierd caracterul lor greoi Nu am sfătui prea mult pe cititori să studieze amănunţit, cu creionul în mînă, exemplul concret împrumutat de la Pacaud, cu diferitele sale variante Intervenţia factorului aleator în duratele operaţiilor presupuse determinate iniţial se reduce în primele două capitole la raţio- narea pe timpi medii şi pe dispersii, care se adună, admiţîndu-se independenţa repartiţiilor şi existenţa nedeplasării în estimaţiile succesive Capitolul al treilea generalizează metoda PERT introducînd un model stochastic, care se aplică, în particular, cazurilor în care nu există nici o experienţă prealabilă, permiţînd stabilirea de grafuri perfect definite Capitolul al patrulea, tratînd optimizarea asociată duratelor şi costurilor, este cel mai dezvoltat şi cel mai pasionant în expunerea algoritmului lui Fulkerson se întîlneşte problema mai generală a programării liniare parametrice şi conceptul universal de dualitate Dacă aspectul teoretic poate să pară puţin, greoi, trebuie ca cititorul, după o primă lectură, să încerce, înainte de a reveni, să urmărească pas cu pas exemplul din 4 9, care ocupă un loc important în lucrare şi constituie unul din elementele cheie Cea mai mare parte a capitolului al patrulea tratează, desigur, optimizarea unui program cu durate fixate ale operaţiilor Doar în ultimele pagini este schiţat un caz aleator, în care timpul operator mediu şi abaterea sa tip variază liniar în funcţie de cost, ceea ce constituie tot ce se poate face în această etapă mai rafinat şi mai aproape de o realitate intrinsec complexă Cazurile mai speciale fac obiectul comunicărilor foarte recente ale specialiştilor, la toate congresele şi în toate revistele periodice de cercetare operaţională Această lucrare, aşa cum se prezintă, poate fi citită şi asimilată de cercul larg al celor interesaţi în aceste probleme Sperăm ca lucrarea să ajute la o mai bună cunoaştere a importanţei şi a supleţei metodei PERT, a cărei notorietate este astăzi foarte mare Astfel, în Statele Unite administraţia ia în considerare numai ofertele pentru executarea cărora candidatul prezintă un program de ordo- nanţare, conform acestei metode în încheiere, voi cita cîteva cuvinte ale directorului general al unei mari societăţi din Spania: ,,De cînd am aplicat recent metoda PERT pentru construirea unei clădiri industriale în uzina noastră, nu mai există nici un proiect nou care să nu facă obiectul stabilirii unui program PERT, în înţelegere strînsă între proiectanţi, inginerii constructori şi cercetătorii operaţionali Dacă am fi cunoscut-o şi aplicat-o mai dinainte, în mod sigur am fi evitat multe ezitări, manevre false, reproşuri între părţile interesate şi am fi cîştigat mult în timp, în bani şi în coeziunea inginerilor noştri“ E VENTURA Director al Societăţii de studii aplicate in cercetarea operaţională PREFAŢĂ Printre metodele elaborate în ultimii ani cu scopul punerii în lucru şi al supravegherii executării programelor de producţie sau de cercetare, metoda drumului critic ocupă un loc de frunte, prin simplitatea şi eficacitatea sa precum şi prin bogăţia variantelor care se dezvoltă neîncetat Fie că este vorba de metoda americană PERT, de metoda franceză a potenţialelor sau de diverse variante propuse cam peste tot, baza matematică a acestora este constituită de teoria grafurilor, ea însăşi o parte importantă a teoriei mulţimilor Este necesar ca inginerii, organizatorii, planificatorii şi, în general, toţi cei care au responsabilităţi în elaborarea şi controlul programelor, să cunoască posibilităţile oferite de metoda drumului critic Trebuie să arătăm că această metodă nu este universală şi că multe probleme nu pot fi rezolvate cu ajutorul ei; pe măsură ce se descoperă şi se introduc noi metode de ordonanţare şi de planificare, se descoperă că apar noi probleme care necesită algoritmi diferiţi sau mai perfecţionaţi Structurile lumii moderne se complică din zi în zi şi metodele ştiinţifice devin indispensabile pentru orientarea în acest gen de dificultăţi combinatorii, în cazul activităţilor mici, mijlocii sau mari în unele ţări metoda PERT este utilizată chiar în întreprinderile mici; cît despre marile întreprinderi, acestea au servicii speciale cu sarcina de a aplica adesea cu ajutorul unor calculatoare electronice puternice, procedeele bazate pe această metodă Sînt numai cîţiva ani de cînd sistemul PERT (Program Evaluation Research Task) a fost utilizat pentru realizarea programului de cercetare şi de construire a cunoscutelor rachete Polaris; drumul parcurs de atunci este considerabil şi nu mai există un număr al unei reviste de cercetare operaţională sau de gestiune ştiinţifică care să nu cuprindă cel puţin un articol asupra acestui subiect, în care să nu fie propusă o nouă dezvoltare Propunem cititorilor acestei lucrări să continuie studierea problemelor de ordonanţare consultînd bibliografia existentă foarte bogată, atît asupra aspectelor teoretice, cit şi asupra celor practice, care îi interesează Din aceste eforturi vor rezulta importante economii de timp şi de bani Cadrul în care se plasează activităţile omeneşti ale timpului nostru nu este numai spaţial sau numai temporal; el este alcătuit din coordonate de spaţiu şi de timp, unde totul este mişcare şi adaptare Matematicienii şi inginerii ne dau posibilităţile navigării în acest cadru Să nu le îndepărtăm prin vanitatea intuiţiei noastre, intuiţie preţioasă, dar pe care trebuie să o ajutăm şi să o amplificăm prin procedee de investigaţie convenabile A KAUFMANN şi G DESBAZEILLE GRAFURISIORDONANTARE 5 5 1 1 INTRODUCERE Fără a face abuz de aspecte teoretice relative la cercetarea relaţiilor ordinale ale unui graf, acesta nefiind altceva, după cum vom vedea, decît o mulţime de elemente între care există legături orientate, credem că este potrivit să amintim anumite noţiuni din teoria grafurilor şi să arătăm cum se utilizează diverse procedee care permit realizarea unor relaţii ordinale Ne vom limita la ceea ce este esenţial pentru un studiu al procedeelor de ordonanţare care intervin în metoda drumului critic Vom defini mai întîi noţiunile de vîrf, arc, drum, circuit etc , apoi vom introduce noţiunea de ordine Scopul acestui studiu preliminar este de a ne învăţa să căutăm în ce ordine pot fi realizate operaţii despre care ştim numai că unele trebuie să le preceadă pe altele 1 2 NOŢIUNEA DE GRAF Să considerăm o mulţime de elemente în număr finit, pe care le vom nota cu literele A1; X2, , Xn, mulţimea ) însăşi notînd-o prin litera X, şi vom scrie simbolic : X = {An X2, ,Xn} (1 1) Dacă fiecărui element din mulţimea X îi corespund zero, unul sau mai multe elemente din X, se va forma, în sensul teoriei mulţimilor, un graf Numind T legea care reprezintă corespondenţa considerată, vom utiliza pentru a reprezenta un graf simbolul G = (X, f) (1 2) Fig 1 2 Fig 1 2 0 1 1 1 0 0 1 1 1 0 0 0 0 1 0 1 1 1 0 0 1 1 1 0 0 0 1 0 1 0 0 0 0 0 0 0 Tabelul 1 1 A B C D E F Să considerăm un exemplu pentru a ilustra această definiţie, «are ne va părea întrucîtva puţin mai abstractă Fie o mulţime de şase elemente şi să considerăm legea următoare de corespondenţă definită prin simbolul T : T2 = {B, C, D}, rB = {A, B, C}, TE, = [B, D, E, F} VD = {C, D, E}, TE = {C, E}, TF =0, unde simbolul 0 reprezintă mulţimea vidă care joacă în teoria mulţimilor rolul lui zero din numeraţie Un graf se poate reprezenta cu ajutorul unui desen ca acela ■din fig 1 1, care se numeşte reprezentarea sagitală a grafului; el se poate reprezenta de asemenea ca în fig 1 2, sau cu ajutorul unui tabel — matrice (tabelul 1 1) în teoria mulţimilor se numeşte graf orice aplicaţie a unei mulţimi în ea însăşi; simbolul T reprezintă această aplicaţie ) Vîrf Un element al unei mulţimi care formează un graf se numeşte -vîrf Pe fig 1 1 elementele A, B, C, D, E şi F sînt vîrfurile grafului Un vîrf este numit de unii autori şi punct Arc O pereche orientată (X;, X,) de două vîrfuri şi X,- se numeşte arc Astfel, în fig 1 1: (A, B), (B, A), (A, (7), (E, E), (C, F) sînt arce Se va nota prin U mulţimea arcelor; în exemplul dat avem : U = {(A, B), (A, C), (A, D), (B, A), (B, B), (B, C), (C, B), (C, JD), (C, E), (C, F), (1 5) (D, C), (D,D), (D, E), (E, C), (E, E)} Un graf poate fi notat prin u = (X, r) sau prin G = (X, U) Drum, Este o succesiune de arce adiacente care permit să se treacă de la un vîrf la altul; astfel (A, B, C, F), (A, C, D, E, C, F), (A, B, A), (A, D), (A, B, B, C) sînt drumuri (fig 1 1) Circuit Este un drum al cărui vîrf iniţial coincide cu vîrful final Astfel (A,C, B, A), (A,C,B,A,C,B,A), (A,B,A), (B, B) sînt circuite (fig 1 1) Buclă Este un arc pentru care originea şi extremitatea coincid; (B, B)f (D, D) şi (E, E) sînt bucle (fig 1 1) Lungimea unui drum sau a unui circuit Este numărul arcelor unui drum sau ale unui circuit Astfel, lungimea drumului (A, B, C,F) este 3, aceea a circuitului (A, D, Ey C, B, A ) este 5 Muchie Cînd două vîrfuri distincte sînt unite prin cel puţin un arc în- tr-un sens şi în altul, se spune că între ele există o muchie Astfel, noţiunea de muchie este o noţiune neorientată a grafului De exemplu, în fig 1 1 : între A şi D există o muchie; între A şi B există o muchie; nu există muchie între B şi D O muchie se va nota (A, D) sau (A, B); mulţimea muchiilor se va nota prin U Astfel, în exemplul nostru (fig 1 1) : U = ((Â7B), (Ă^), (Ă^D), (B^C), (CJ)}, (CTE), (CŢF), (&E)} (1 6) Lanţ Este o succesiune de muchii adiacente, adică o succesiune, de arce adiacente unde nu ţinem seama de orientarea arcelor Astfel, pe fig 1 1 : (B, C, A, D, E) nu este un drum (din cauza arcului AC), însă este un lanţ Orice drum este evident şi un lanţ; însă, reciproc, după cum am văzut, un lanţ nu este totdeauna un drum Graf conex Un graf este conex dacă între orice pereche de vîrfuri ale grafului există cel puţin un lanţ Graful din fig 1 1 este conex Un graf este tare conex dacă între oricare pereche de vîrfuri există cel puţin un drum Graful din fig 1 1 nu este tare conex (nu există drumuri de la F către celelalte vîrfuri) Graf antisimetric Să considerăm mulţimea U a arcelor grafului ; dacă pentru orice arc (X4, X,) aparţinînd lui U arcul (X,,Xj) în sens invers nu aparţine lui U, se spune că graful este antisimetric Rezultă din această definiţie că un graf antisimetric nu are bucle De exemplu, graful din fig 1 3 este antisimetric Graf fără circuite în prezenta lucrare ne vom o- cupa în special de grafuri fără circuite, adică neconţinînd nici un circuit, aşa cum este cazul grafului din fig 1 3 Este evident că orice graf fără circuite este antisimetric, însă reciproca nu este adevărată amintit numai ceea ce este strict necesar pentru o prezentare logică a conceptelor pe care le vom introduce în metoda drumului critic, obiectul acestei cărţi amintit numai ceea ce este strict necesar pentru o prezentare logică a conceptelor pe care le vom introduce în metoda drumului critic, obiectul acestei cărţi Cititorul care cunoaşte teoria grafurilor va fi uimit de faptul că ne limităm la aşa de puţine lucruri dintr-o teorie atît de amplă; am Tabelul 1 2 Matricea corespunzătoare grafului reprezentat in fig 1 3 Zerourile nu au fost trecute în această matrice RELAŢIA DE ORDINE STRICTĂ ÎNTR-UN GRAF CONEX FĂRĂ CIRCUITE Să considerăm relaţia următoare între două vîrfuri distincte X4 şi Xj ale unui graf conex fără circuite : există un drum de lungime nenulă de la Xt spre X,, pe care o vom nota X4 —► a 1 1 0 X X b 1 1 2 1 0 X c 0 X X X d 1 1 0 X X e 1 1 1 0 X f 0 X X X g 1 1 2 2 2 0 h 0 X X X i 1 1 2 1 0 X i 1 1 0 X X k 1 1 1 0 X l 1 1 1 0 X m 1 1 2 2 2 0 1 5 GĂSIREA UNEI ORDONANŢĂRI ) ORDONANŢARE DE VÎRFURI Vom spune că intr-un graf conex fără circuite, am realizat o ordonanţare a vîrfurilor dacă am construit un şir al tuturor vîrfurilor grafului care satisfac toate condiţiile de anterioritate definite de arcele grafului Un astfel de şir nu este deci altceva, decît o permutare a vîrfurilor satisfăcînd toate condiţiile de anterioritate Pentru a construi o ordonanţare a vîrfurilor este suficient să alegem vîrfurile respectînd ordinea nivelelor definite în 1 4 Exemple în figurile 1 3 şi 1 4 : EJIKACDHFBG, IJEAKHCFDGB sînt ordo- nanţări ale vîrfurilor grafului în figurile 1 7, 1 15 şi 1 16 : HEIAKDJCBGFML (fig 1 15) Şi HE1ADKJBCFL (fig 1 16) sînt ordonanţări ale vîrfurilor grafului din fig 1 7 Vom observa că procedeele pe care le-am descris pentru a găsi nivelele nu ne permit, aşa cum sînt date, de a enumera toate permutările care satisfac condiţiile de anterioritate date de graf Pentru ceea ce ne preocupă în lucrarea aceasta, o astfel de enumerare nu prezintă interes, ea puţind prezenta totuşi un interes evident în alte probleme sau alte metode ORDONANŢAREA ARCELOR Considerînd un graf conex fără circuite, vom spune că am realizat o ordonanţare a arcelor dacă am construit un şir din toate arcele, satisfăcînd toate condiţiile de anterioritate ale arcelor, unele în raport cu celelalte Un astfel de şir nu este altceva decît o permutare a arcelor care satisfac toate condiţiile de anterioritate ale arcelor Pentru a construi o astfel de ordonanţare este suficient de a alege arcele,, respectînd ordinea nivelelor arcelor definite în paragraful precedent Exemplu (fig 1 17 şi fig 1 21) chf dj akbleimg şi hcfadjbkliegm sînt ordonanţări ale arcelor grafului din fig 1 17 ORDONANŢAREA COMPLETĂ A UNUI GRAF CONEX FĂRĂ CIRCUITE Se va spune că s-a realizat o ordonanţare completă a unui astfel de graf dacă, după ce s-au desemnat simbolurile pentru reprezentarea vîrfurilor şi simbolurile pentru reprezentarea arcelor, s-a realizat pe de o parte o ordonanţare a vîrfurilor şi pe de altă parte o ordonanţare a arcelor Exemplu Să considerăm graful fără circuite din fig 1 22 Am arătat cum se determină nivelele; o ordonanţare a vîrfurilor va fi, de exemplu: 1, 2, 3, 6, 4, 8, 5, 7, 11, 9, 10, 12 Tabelul 1 6 1 2 3 4 5 6 7 8 9 10 11 12 —► Fo “► —> V, -► F* —► v -> ra 1 1 1 1 3 3 3 3 3 3 2 1 0 2 1 1 1 3 3 3 3 2 2 1 0 X 3 1 1 1 1 4 4 4 3 3 2 0 X X 4 1 1 1 3 3 3 2 1 0 X X X 5 1 1 1 1 0 X X X X X e 1 1 1 3 3 3 2 1 0 X X X 7 1 1 1 1 0 X X X X X 8 1 1 1 1 4 4 3 1 0 X X X X 9 1 1 1 0 X X X X X X 10 1 1 0 X X X X X X X 11 1 1 2 1 0 X X X X X X 12 0 X X X X X X X X 12 10 9 5 8 4 3 2 1 11 6 6 Rămîne de găsit o ordonanţare a arcelor Operînd ca mai înainte, la 1 4, cititorul va găsi uşor cum a fost obţinută fig 1 24 prin eliminarea progresivă a arcelor fără antecedenţi O ordonanţare a arcelor va fi, de exemplu : aebcfghjndvmilopqksutrwzyx în anumite cazuri ne va interesa o descompunere a nivelelor în subnivele, unde sînt grupate arcele avînd acelaşi vîrf ca origine (fig 1 25) Algoritmul dat în 1 4 pentru determinarea unei ordo- nanţări a arcelor va fi uşor modificat în scopul obţinerii de subnivele ale arcelor avînd o origine comună UTILITATEA STUDIULUI UNEI ORDONANŢĂRI COMPLETE După cum se va vedea mai departe, va trebui să ne preocupăm de programe cuprinzînd operaţii foarte diferite şi adesea date şi specificate „de-a valma” fără a se şti bine care sînt acelea care trebuie să se plaseze înaintea altora, cunoscîndu-se numai anumite antecedente obligatorii De fapt este vorba, în primul rînd, de a construi graful reprezentînd relaţiile cunoscute din program, adică antecedentele ; acest lucru fiind făcut, se întreprind apoi alte cercetări importante Scopul acestui prim capitol a fost tocmai de a arăta cum se poate construi graful reprezentativ şi cum se pot numerota convenabil vîrfurile şi arcele Anumite programe formează grafuri avînd mai multe mii de vîrfuri şi mai multe mii de arce, iar tratarea lor la calculatorul electronic implică mai întîi cunoaşterea unei oidonanţări complete Este un principiu de organizare după care ne conducem în astfel de consideraţii : tot ceea ce poate fi ordonat trebuie să fie ordonat Vom examina acum grafuri în care sînt afectate arcelor valori numerice (intervale de timp, termene etc ) STABILIREA UNUI PROGRAM DE CERCETARE SAU DE PRODUCŢIE 2 1 INTRODUCERE Vom arăta în acest capitol cum se stabileşte graful unui program de cercetare, de dezvoltare, sau de producţie Apoi, prin calculul drumului critic vom deduce timpul de realizare a ansamblului de lucrări, intervalele de variaţie şi marginile operatorii, acestea în cazul cînd operaţiile au durate bine determinate Vom avea în vedere, mai departe, cazul în care variabilele reprezentînd duratele operaţiilor sînt aleatoare, ceea ce constituie cazul considerat în metoda PERT propriu-zisă Vom da diverşi algoritmi permiţînd calcularea drumului sau drumurilor critice şi diverselor abateri cărora li se ataşează o mare importanţă în sfîrşit, vom vorbi de tratarea acestor programe cu ajutorul calculului electronic 2 2 REPREZENTAREA UNUI PROGRAM PRINTR-UN GRAF 2 2 1 GENERALITĂŢI în ceea ce va urma, vom numi program o mulţime de operaţii care conduc la realizarea unui obiectiv (proiect, mare ansamblu «tc ), astfel încît se cunosc pentru fiecare dintre operaţii durata sa (determinată sau aleatoare) şi relaţiile de ordine privitoare la aceasta (antecedente obligatorii) Un program poate fi reprezentat printr-un graf, numit graf de program sau graf de ordonanţare, ale cărui arce reprezintă operaţiile ; duratele operaţiilor sau timpii operatori fiind presupuşi cunos- cuţi, se va ataşa arcelor grafului o valoare nenegativă ) reprezen- tînd fiecare din aceşti timpi operatori Vîrfurile grafului, numite evenimente şi notate Ef, pot fi interpretate ca indicînd realizarea obiectivelor parţiale (montarea unui subansamblu, de exemplu) Dacă fineţea analizei (operaţii din ce în ce mai simple şi cu durate din ce în ce mai scurte) este foarte mare, evenimentele considerate pierd mult din semnificaţia lor practică, îndeosebi unele dintre ele, care indică marile etape ale proiectului Durata operaţiei Po, care este una din operaţiile contribuind la realizarea evenimentului E}, va fi notată tiS Figura 2 1 reprezintă un program Stabilirea grafului necesită pentru fiecare operaţie cunoaşterea operaţiilor care o preced imediat Graful traduce astfel relaţiile de ordine existente în mulţimea operaţiilor; el constituie o reprezentare a ceea ce unii numesc ordonanţare ) Se va observa că unei aceleiaşi ordonanţări a operaţiilor, deci aceluiaşi graf, îi pot corespunde mai multe programe diferite, dacă timpii operatori sînt diferiţi Se va observa că graful unui program este fără circuite, altfel o operaţie ar putea fi făcută să urmeze ei însăşi Primul lucru de îndeplinit în vederea aplicării metodei drumului critic va fi, astfel, întocmirea unei liste de operaţii care va fi subdi- vizată, mai mult sau mai puţin, în raport cu precizia dorită Ne vom strădui să alegem operaţii cit mai simple, pentru care timpii operatori să fie de valori comparabile Apoi se vor determina antecedentele fiecăreia dintre operaţii Această parte numită analiză de lucru este cea mai importantă şi adesea cea mai lungă Ea necesită colaborarea strînsă a întregului personal (ingineri, contramaiştri, muncitori) interesat în realizarea pe care vrem să o întreprindem şi care se dovedeşte foarte utilă ca atare Pentru a concretiza, vom examina un exemplu, acela al construcţiei unei vile Pentru a construi o vilă este necesar un anumit număr de operaţii pe care le enumerăm într-o anumită ordine : elaborarea planurilor; obţinerea autorizaţiei de construire ; semnarea contractului; transportarea materialului la şantier (agregate, piatră de construcţie, cărămizi ); executarea racordurilor electrice şi de apă pentru şantier; executarea fundaţiilor; executarea lucrărilor de zidărie portantă; executarea şarpantelor; executarea acoperişului; executarea pardoselilor; executarea pereţilor despărţitori din cărămizi uşoare; executarea instalaţiilor de încălzire centrală, gaze ; executarea instalaţiilor electrice interioare ; montarea tocurilor; tencuirea pereţilor; executarea a două şemineuri în piatră; montarea obiectelor sanitare şi a bucătăriei; executarea lucrărilor în ipsos ; montarea prizelor şi a întreruptoarelor electrice ; montarea uşilor şi a ferestrelor; executarea zugrăvelilor şi a tapetelor; montarea centralei termice şi a radiatoarelor; transportarea pămîntului pentru grădină; montarea contoarelor şi executarea branşamentelor pentru energie electrică, gaze, telefon, apă; curăţirea şantierului; recepţia provizorie Să numim E1 evenimentul iniţial sau începutul lucrărilor şi să numim E2 evenimentul cînd contractul poate fi semnat între E± şi E2 se găsesc două operaţii: elaborarea planurilor şi obţinerea autorizaţiei de construire Dacă regulamentul administrativ prevede că nu se poate acorda o autorizaţie de construire fără a avea planurile, cum se întîmplă în general, vom aşeza între Et şi E2 un eveniment E3: planurile sînt terminate După evenimentul E2, se vor putea întreprinde simultan operaţiile e şi c, cu timpii operatori t24 şi Z25; într-adevăr, racordurile electrice şi de apă nu sînt în general făcute de antreprenor Evenimentul E5 constă în semnarea contractului; se poate atunci întreprinde operaţia d, de timp operator Z54 Operaţia/, de timp operator 7 unde (i, j) reprezintă arcul mergînd din vîrful i la vîrful j şi Uf este submulţimea arcelor care se întîlnesc în j şi sînt orientate către j Fig 2 28 Fig 2 2J Fig 2 28 Fig 2 2J MAX înseamnă că se ia maximul pentru toate arcele (i, j) apariţinînd lui Uj' De asemenea: : t; = MIN (t) - tj, i = 2, , n - 1, (2 6) t„ = t„ — data de realizare a proiectului total, (2 7) t‘i = ) g V fiind mulţimea drumurilor mergînd de la j către n, - o, (2 17) formulele (2 9) — (2 12) determină deci intervalele de fluctuaţie şi marginile Algoritmii lui Ford şi Bellman-Kalaba vor fi de asemnea folosiţi în această căutare a drumului critic şi a datei tn de realizare a programului Iată, cu titlu de informaţie, aceşti doi algoritmi ALGORITMUL LUI FORD 1° Se notează fiecare vîrf Et cu tt; se începe prin a lua t± = 0 şi ti = 0 pentru i^= 1 2° Se caută un arc (/, j) astfel încît t, — tj; se va observa că t) > 0 dacă j =f= 1 Se continuă astfel pînă cînd nici un arc nu mai permite să-l mărim pe pentru orice i 1 3° Există un vîrf EPl astfel încît tm tI>l - tP1 n, căci tn a crescut monoton în cursul procedeului şi EPl este ultimul vîrf folosit pentru a-l mări pe t„ La fel fie EPî astfel încît tPl tPî = tPlPi Şirul tn, tp„ tPl, fiind strict descrescător, se va ajunge la un anumit moment la situaţia EPlc + r = Ev Atunci tn va fi valoarea drumului cel mai lung de la Ek la En (data de realizare a programului) şi acest drum va fi p = [Uj E„k, En v , EPi, E„J ALGORITMUL BELLMAN-KALABA Este o variantă a celui precedent prezentat cu ajutorul programării dinamice El se bazează pe proprietatea evidentă următoare : orice drum maximal cuprinzînd cel mult r arce este format din drumuri parţiale conţinînd cel mult k arce (k C r) şi care sînt de asemenea maximale *) Pentru orice arc (i, j) e II se pune = — oo; pe de altă parte pentru orice i se pune tti = 0 Se propune găsirea unui drum l1 — [ -®1, • • • > -®nl, (2-18) astfel incit ^li, + + tiii, + • • • + n (2-19) să fie maximal E6te suficient să rezolvăm sistemul de ecuaţii: Vi = MAX + v,), i = 1,2, , n — 1, (2 20) = 0,’ (2 21) de la 1 la n originea fiind vîrful 1 şi terminaţia vîrful n Se operează după cum urmează, punînd mai întîi: ®T = «in, i = l,2, , n-l, (2 22) fl'»1’ = 0 (2 23) Se calculează apoi: = MAX (ti, + t,',1’), i = 1, 2, , n — 1, ' j = l,2, ,n, (2 24) v™ = 0 ; (2 25) de la 1 la n originea fiind vîrful 1 şi terminaţia vîrful n Se operează după cum urmează, punînd mai întîi: ®T = «in, i = l,2, , n-l, (2 22) fl'»1’ = 0 (2 23) Se calculează apoi: = MAX (ti, + t,',1’), i = 1, 2, , n — 1, ' j = l,2, ,n, (2 24) v™ = 0 ; (2 25)unde Vi (i = 1,2, , n — 1) reprezintă valoarea drumului optimal de la vîrful i la vîrful final Cele n vîrfuri sînt presupuse numerotate ♦) Această proprietate este de asemenea adevărată dacă se consideră drumurile minimale şi drumurile parţiale minimale apoi succesiv : r'î’=MAX i = l,2, ? j = 1, 2, n, (2 26) v = 0, (2 27) pentru toate valorile lui k = 1,2,3,4,- ; ne oprim cînd i = 1, 2, , n (2 28) Atunci v(î’ reprezintă valoarea drumului optimal între vîrfurile Et şi E„ Se demonstrează uşor că dacă graful are n vîrfuri, sînt suficiente n—2 iteraţii pentru a atinge optimul Este important de subliniat că algoritmul Bellman-Kalaba stabilit la început pentru căutarea unui drum cu valoarea minimală se aplică, aici, la căutarea unui drum cu valoare maximală, căci graful unui program nu conţine nici circuite, nici bucle, în caz contrar, lungimea, numărul arcelor conţinute de un drum putînd să nu fie finită, iar valoarea acestui drum nu ar fi mărginită EXEMPLU DE APLICARE A ALGORITMULUI BELLMAN-KALABA Fie graful din fig 2 30, unde numerele înscrise pe arce sînt timpii operatori Fig 2 30 Tabelul 2 1 1 2 3 4 6 6 7 8 9 10 11 12 1 0 3 4 11 — 00 8 — 00 — 00 — oo — 00 — 00 — 00 2 — 00 0 — 00 7 9 — oo 7 — 00 — 00 00 — — 00 — 00 3 — 00 — 00 0 9 — 00 2 4 — 00 — 00 — 00 — 00 — 00 4 — 00 — 00 — 00 0 2 ~ 00 — 00 3 3 — 00 — 00 — 00 6 — 00 — 00 — 00 — 00 0 — 00 — oo 5 — 00 — 00 — 00 — 00 6 — 00 — 00 — 00 — 00 — 00 0 5 — 00 8 9 — 00 — 00 7 — 00 — 00 — 00 — 00 — oo — 00 0 4 — 00 8 5 — 00 8 — 00 — 00 — 00 — 00 — 00 — 00 — oo 0 — 00 — oo 11 — 00 9 — 00 — 00 — 00 — 00 — 00 — 00 — oo — 00 0 4 — 00 5 10 — 00 — 00 — 00 — 00 — oo — oo — 00 — 00 — 00 0 4 3 11 — 00 — 00 — 00 — 00 — 00 — 00 ~ 00 — 00 — oo — 00 0 2 12 — 00 — 00 — 00 — 00 — 00 — 00 — 00 — 00 — 00 — 00 — 00 0 Pentru a calcula valorile ale drumurilor optimale de la E1 la E12, s-a întocmit tabelul 2 1,, care reprezintă matricea valorilor asociate grafului, apoi s-au calculat , prin formulele (2 22) — (2 28) Se obţine tabelul 2 2 Tabelul 2 2 \ i k \ 1 2 8 4 5 6 7 8 9 10 11 12 1 — 00 — 00 — 00 — 00 — 00 — 00 — 00 — 00 5 3 0 2 — 00 — 00 — 00 8 — 00 13 11 |13| 7 Î6Î 3 21 18 17 16 |18| 16 |17| 13 |îoj 6 2 0 4 27 l27l 25 |20| 18 22 17 13 10 6 2 0 5 31 27 20 18 22 17 13 10 6 2 0 6 I33 27 29 20 18 22 17 13 10 6 2 0 7 33 27 29 20 18 22 17 13 10 6 2 0 Să dăm detaliat o parte din calcul; să vedem, de exemplu, cum s-a obţinut = 16 = MAX + ® ] = = max [Z61 + v'i}, t62 + v{?, tn + v{?, t6i + v{l', t65 + v{l}, r66+®(02), ^67 + ^68 + ^69 + 'D 2)]}, (^a) — max I MAX [v2 (r2) -f- ®23 (r2, x3) ] , l (2 29) MAX (ajj) -ţ- r13 (®j, r3)] , $93 (r0, x3) | ' P î } | L l Stadiu Instalaţie / / > y z / / X /Pstratulde protecţie in uzine ţ/stropiri' - Grinda^ ; Studiu ' Instalaţii şantier Amplasare D o \Fufida!ti debelon Grinzi de betan Pordoreah cudaie/ Perete de umplutura Tencuială t Studiu \ J' Construcţia in atelier Transport \c ’k v v A / — — - '• ” y ' Z / 7 z / Studiu \ 1 Executarea izolaţiei A i \ / X Z | Studiu \ Montaj / I ■ 9-—>' P î } | L l Stadiu Instalaţie / / > y z / / X /Pstratulde protecţie in uzine ţ/stropiri' EundoSe — 'G EundoSe — 'G Timpii operatori sînt indicaţi în tabelul 2 4 (în zile) Putem trasa, atunci, graful programului (fig 2 35) Tabelul 2 4 A - 10 H - 20 B - 35 B - 3 I - 4 S - 40 C - 15 J - 18 T - 2 D - 10 K - 60 U - 20 E - 15 L — 20 V - 20 F - 5 M - 5 IV - 35 G - 30 N - 25 X - 55 Au trebuit introduse două operaţii virtuale care sînt punctate pe figură, Numerele notate în vîrfurile grafului sînt datele aşteptate ale evenimentelor corespunzătoare Vîrful terminal poartă data 250 Ansamblul lucrărilor va fi deci terminat în 250 de zile Drumul critic este figurat prin trăsături groase Dar anumite operaţii pot fi descompuse ; altfel spus, pot fi întreprinse acele operaţii care urmează înainte de terminarea completă a precedentelor Aceasta se va produce, dacă clădirea se construieşte travee cu travee Diferite programe pot fi atunci examinate, fiecare dintre ele corespunzînd unei organizări diferite a şantierului 1° Programul obţinut construind cu prioritate prima travee Operaţiile descompuse sînt A, C, K, L, M, N, P, S, U, V Desemnarea operaţiilor şi timpilor operatori sînt indicate în tabelul 2 5 Operaţiile Operaţiile A2 Studiu pierderi, alegere aeroterme, traseu, deviz greutate suspendată A2 Studiu execuţie a instalaţiilor de încălzire B Studiu acoperiş Cj Notăm calcule, deviz greutate şarpantă C2 Studiu execuţie şarpantă D Studiu beton armat E Terasamente generale F Amplasare, instalaţii şantier G Terasamente fundaţii H Beton fundaţii 1 Terasamente grinzi J Beton grinzi Kj Fabricarea şarpantei primei travee lemn K2 Fabricarea celorlalte patru travee Lj Strat protecţie al primei travee jL2 Strat protecţie al celorlalte pentru travee A2 Studiu pierderi, alegere aeroterme, traseu, deviz greutate suspendată A2 Studiu execuţie a instalaţiilor de încălzire B Studiu acoperiş Cj Notăm calcule, deviz greutate şarpantă C2 Studiu execuţie şarpantă D Studiu beton armat E Terasamente generale F Amplasare, instalaţii şantier G Terasamente fundaţii H Beton fundaţii 1 Terasamente grinzi J Beton grinzi Kj Fabricarea şarpantei primei travee lemn K2 Fabricarea celorlalte patru travee Lj Strat protecţie al primei travee jL2 Strat protecţie al celorlalte pentru travee Tabelul 2 5 Tabelul 2 5 Timpii operatori Timpii operatori 3 3 10 A 15 C 10 A 15 C 10 15 5 30 20 4 10 15 5 30 20 4 18 18 60 60Transport prima travee M2 Transport celelalte patru travee Montaj ridicare prima travee Nt Montaj ridicarea celorlalte patru travee Montaj instalaţii de Incălzirre prima travee Montaj celelalte patru travee Sj Zugrăveli prima travee S2 Zugrăveli celelalte patru travee T Studiu instalaţii electrice Uj Izolaţii acoperiş prima travee U2 Izolaţii acoperiş celelalte patru travee Instalaţii electrice prima travee V2 Instalaţii electrice celelalte patru travee W Pardoseală cu dale 8 32 2 4 16 8 32 2 4 16 S S U U 35 55 35 55X Tencuirea pereţilor de umplutură Acum se poate trasa graful acestui nou program şi calcula datele aşteptate ale diverselor evenimente ale programului (fig 2 36) 8e vede că data evenimentului final este acum 210 S-au cîştigat deci 40 zile descompunînd anumite operaţii Ca şi mai înainte, operaţiile punctate sînt operaţii virtuale, a căror motivare a fost indicată în 2 2 2° Programul obţinut construind cu prioritate două travee Noii timpi operatori sînt : Tabelul 2 6 K1 - 24 N, - 10 17,-8 K„ - 36 N2 - 15 Ui - 12 — 8 - 14 V,- 8 L2 - 12 R2 - 21 V2 - 12 JW, - 2 S, - 16 M2 - 3 S2 - 24 Ceilalţi timpi operatori rămîn neschimbaţi Se obţine graful din fig 2 37 Data evenimentului final este 177 S-a redus deci din nou durata ansamblului lucrărilor 3° Programul obţinut construind cu prioritate trei travee Noii timpi operatori sînt: Tabelul 2 7 «i- 36 N, - - 15 ut - - 12 *2 - 24 A» - - 10 u2 - - 8 il - 12 - - 21 V!- 12 l2 - 8 «2 - - 14 v,- - S - 3 Si - - 24 m2 - 2 s2 - - 16 Ceilalţi timpi operatori rămîn neschimbaţi Se obţine graful din fig 2 38 Data evenimentului final este 165 Operaţiile critice sînt redate cu linii groase Am prezentat trei programe corespunzînd la trei organizări posibile ale şantierului, însă mai există şi altele Interesul acestui exemplu era de a da o privire generală asupra a ceea ce s-a convenit să se numească probleme de ordonanţare şi de a arăta cum metoda drumului critic permite să alegem între mai multe programe ale aceluiaşi ansamblu de lucrări în exemplul nostru, vom alege programul trei care dă un termen total de execuţie de 165 zile 2 C 2 AL DOILEA EXEMPLU Acest nou exemplu ne-a fost comunicat de B Roy şi a fost extras dintr-un studiu făcut de Rosinski ), exemplul ilustrîndu-ne importanţa metodei potenţialelor în problemele de ordonanţare, la fel cu uşurinţa sa de aplicare Problema pusă este aceea a realizării ciclului zilnic pentru o nouă unitate de tunele-cofraje concepută în vederea realizării fundaţiilor unui microraion de locuinţe Să precizăm, mai întîi, că operaţiile de grutaj necesită în total aproape 10 ore pe ciclu şi că între ele trebuie intercalate anumite lucrări manuale (9 ore de lucru pentru un muncitor în medie); se concepe că realizarea unui plan executabil în mai puţin de 12 ore ar fi delicată Figura 2 39 reprezintă graful asociat acestei probleme Datele (durate, lungimi de arce) sînt relative la cadenţele de lucru minimale (semnul -f- indică duratele şi arcele „comprimabile”) Se vor nota printr-o încadrare dublă sarcinile recurgînd la macara Se va observa de asemenea prezenţa restricţiilor de formă particulară între sarcinile 104 şi 106 exprimînd că acestea sînt legale în execuţia lor (106 trebuie să înceapă cu 25 de minute după 104r nici mai înainte, nici mai tîrziu); este una din posibilităţile formulării prin metoda potenţialelor Lăsăm de o parte din cauza spaţiului descrierea completă a operaţiilor Pe aceste baze, se va putea aplica unul din algoritmii (din 2 5) grafului din fig 2 39 şi se vor obţine astfel datele de început ale primelor sarcini care figurează pe acest graf Se constată imediat contradicţii: de exemplu, sarcinile 23 şi 51 pe de o parte, 6, 9 şi 10 pe de altă parte (şi încă altele) fac apel simultan la macara, ceea ce este imposibil O regulă de arbitraj de felul ,,primul venit — primul servit” permite eliminarea acestor contradicţii şi dă adesea soluţii apropiate de optim, procedînd după cum urmează: calculul se face cu ajutorul algoritmului ales, dar clasificînd momentele' de început Z, în ordine crescătoare, dacă există o disjuncţie între două sarcini i şi j (Z,- ZJ, ea elimină cu ajutorul unei restricţii potenţial suplimentare : unde dj este durata sarcinii i, ea impunînd lui j să urmeze după i Putem urmări procesul pe graful figurii 2 40 care dă prima soluţie : sarcina 51 trebuie să urmeze sarcinii 23 (Z51 = t22, arbitraj asupra numărului sarcinii) şi trebuie să preceadă sarcinii 2 (35+) DATELE DE BAZA p 7 Data de începere Sarcina nr 52 ■ ■ (17+) Duratamaxima 17minute,comprimaiile' II || Sarcini cu grutaj 3U Valoarea iui ti 'data cea mai apropiata I 571 447 Valoarea lui tpdata cea mai tirue — *- —Restricţii disjunctive Drum critic Calculul datelor limită/* sau al datelor „cele mai tîrzii” se face în continuare, fără a ţine seama de disjuncţii, pornind de la graful prevăzut cu arce disjunctive (punctate pe figură) Toate datele nu figurează pe desen pentru a nu-l încărca prea mult Această primă soluţie are deja o durată inferioară lui 12 h (684 minute) Trecerea de la o soluţie la alta mai bună se face prin revizia şi explorarea sistematică: a ipotezelor asupra cadenţelor de lucru, a modului de tratare a restricţiilor disjunctive Pentru a comprima durata totală, trebuie accelerată execuţia sarcinilor drumului critic Eliminarea unei disjuncţii printr-o restricţie potenţial este o decizie arbitrară, şi se poate întotdeauna lua în consideraţie inversarea ordinii astfel impuse Apriori, aceasta conduce la încercarea unui număr foarte mare de combinaţii Pentru a elimina sistematic pe cele mai puţin bune, se poate utiliza regula următoare : parcurgînd arcele drumului critic de la origine, să se inverseze primul arc disjunctiv întîlnit; să se refacă calculul acestui nou program fără a se atinge de celelalte arce disjunctive Acest procedeu se justifică cum urmează : inversarea unui arc nesituat pe drumul critic nu l-ar scurta pe acesta, căci el nu împrumută acest arc; ea poate însă să-l lungească şi arcul devine atunci critic; inversarea unui arc, eliminînd anumite disjuncţii, poate autoriza un program mai scurt; ar conveni deci, în acelaşi timp, să se elimine restricţiile disjunctive devenite inutile sau acelea care ar putea să se opună unei astfel de ameliorări Dar, în acest din urmă caz, arcele despre care este vorba devin critice şi vor fi deci reexaminate într-un stadiu următor Se vede că acest procedeu simplu nu elimină nici o soluţie, bună, dar nu ameliorează programul la fiecare iteraţie El permite totuşi să se cîştige timp şi să obţinem rezultate bune în acest caz Pentru prima soluţie, sarcina 14 — lungă dar puţin importantă — ar întîrzia ansamblul operaţiilor Ar fi suficient să dublăm efectivul corespunzător pentru a cîştiga 44 minute Dar în soluţia 2 arcele disjunctive legînd sarcinile 23,51 şi 2 ar deveni critice; inversarea arcului 51 -> 2 ar da soluţia 3 Practic, această soluţie ar fi cea mai rapidă: arbitrajele disjunctive situate dincolo de sarcina 15 ar fi sensibil optimale, şi rarele sarcini comprimabile situate pe drumul critic ar interveni la momente cînd încărcarea personalului ar fi vizibil foarte mare Dar utilizatorul alege în definitiv o soluţie mai puţin bună: pentru o pierdere de cîteva minute, este preferabil pentru motive de securitate să se execute împreună şi mai devreme sarcinile 22, 25 şi 26 (deplasarea şi aşezarea pasarelelor mobile de circulaţie) Figura 2 41 dă soluţia reţinută: durata totală s-a redus la 639 minute Obiectivul căutat — realizarea ciclului zilnic — ar putea fi deci atins, cu un efectiv inferior celui folosit înainte şi aceasta cu preţul unei investiţii mici faţă de cheltuielile în joc 2 7 OPERAŢII CU DURATA ALEATOARE Cititorul a luat cunoştinţă acum de aspectele principale pe care le oferă utilizarea unui graf pentru reprezentarea unui program de cercetări sau de producţie, în cazul în care duratele operaţiilor sînt cunoscute şi sigure Importanţa unei astfel de aproximări a unui program de cercetări sau de dezvoltare este evidentă : — Planificatorul sau organizatorul au, datorită folosirii grafului, o prezentare clară şi simplă a programului în care apar interacţiunile diferitelor lucrări Stabilirea grafului este deci, în sine, o primă fază foarte importantă, care permite să se precizeze desfăşurarea operaţiilor — Punerea în evidenţă a unuia sau a mai multor drumuri critice permite determinarea operaţiilor care condiţionează desfăşurarea ansamblului programului Aceste aplicaţii critice vor trebui deci, în cursul realizării, să fie supravegheate foarte de aproape şi asupra lor se va îndrepta efortul principal — Operaţiile diferite de cele ale drumului critic apar, din contră, mai puţin rigide şi tolerează anumite întîrzieri relative la datele lor de punere în execuţie (intervale de fluctuaţie, margini libere) In practică, în foarte multe programe în primul rînd la cele de cercetare şi de dezvoltare, duratele ale operaţiilor sînt insuficient cunoscute şi nesigure Se pot prezenta atunci două cazuri: sau operaţiile, despre care este vorba, nu sînt noi şi se cunoaşte pentru fiecare dintre ele într-un mod aproximativ legea de repartiţie a probabilităţilor timpului operator, sau aceste operaţii sînt în întregime noi, sau foarte puţin cunoscute şi nu se cunoaşte legea de repar- Fig 2 41 a ORDONANŢARE OPTIMALA F 18 166 18 166Valoarea lui Valoarea reţinută (daca Ţ t 'L) Drum critic Fig 2 41 q tiţie a fiecărui timp operator Primul caz corespunde, în general, programelor industriale de dezvoltare, al doilea — programelor de cercetare Dacă se cunosc repartiţiile timpilor operatori, nu există dificultăţi pentru determinarea duratei medii (speranţa matematică,) şi a dispersiei fiecărei operaţii Se obţine astfel: *u Şi ®?r Dacă nu se cunosc distribuţiile timpilor operatori, se presupune, în metoda PERT, pentru motive de comoditate în calcule, că aceşti timpi operatori urmează o lege de repartiţie p *) Nu se poate justifica într-un mod riguros alegerea acestei legi, avînd o funcţie de densitate de probabilitate de tip Gauss-Laplace, dar este mărginită la stînga şi la dreapta Figura 2 42 dă curba care reprezintă o lege de probabilitate p Legea p considerată în metoda PERT a fost aleasă astfel încît: E(t) = t = — (J + 4Jf + B), (2 32) 6 ') Vezi în anexă, o relatare privind această lege A, B fiind marginile inferioară şi superioară ale variabilei aleatoare tf M fiind valoarea cea mai probabilă sau moda Specialiştii responsabili cu fiecare operaţie sînt chestionaţi, pu- nîndu-li-se următoarele trei întrebări: La cît estimaţi durata minimală a operaţiei PiS (estimaţie optimistă) î Fie ai} această durată minimală La cît estimaţi durata maximală a operaţiei Ptj (estimaţie pesimistă) î Fie bti această durată maximală Care este durata cea mai probabilă a lui Ptj (moda) î Fie wii3 durata cea mai probabilă Astfel de informaţii fiind subiective se convine să fie strînse cu anumite precauţii în particular, se stabileşte să se realizeze legături între diferitele informaţii provenite de la specialiştii într-o aceeaşi operaţie (de exemplu, un şef de laborator şi adjuncţii săi) care ar trebui eliberate de orice forme de presiuni exterioare, care ar putea deplasa estimaţiile Deoarece nu cunoaştem nimic, vom presupune că aceste informaţii constituie estimaţii nedeplasate ale lui AiJt BtJ, Mi} Fie: = = (2 34) Putem scrie atunci A tii A tii A;, -|- 4- B{j A;, -|- 4- B{j (2 35) (2 35) Notaţia to are semnificaţia că valoarea celui de-al doilea membru al ecuaţiei (2 35) constituie o estimaţie nedeplasată a valorii medii t(j; tif este un estimator al lui t0 Din punct de vedere teoretic, nu avem dreptul să scriem (2 36) (2 36) consideraţie aproximaţiile făcute (asimilarea estimaţiilor subiective cu estimaţiile obiective) şi ipotezele complet empirice admise mai înainte, nu este cu totul greşit de a lua ca estimaţie a dispersiei valoarea 2 = 2 = (2 37) (2 37)Am obţinut aBtfel, cu un anumit număr de concesii la rigoare, esti- maţii ale duratei medii a operaţiei Pj3- şi a dispersiei timpului său operator Vom reţine aceste estimaţii şi vom scrie suprimînd bara şi / Bau accentul circomflex : (2 38) Avînd astfel pentru fiecare operaţie valoarea medie şi «dispersia uf, a timpului operator, fie că sîntem în cazul în care se cunoaşte legea de repartiţie, fie că sîntem în cazul empiric, cînd se posedă numai informaţii apriori, se foloseşte metoda din 2 3 pentru determinarea valorilor medii ale evenimentelor Et şi a drumului critic (sau a drumurilor critice), afectînd fiecărui arc valoarea medie tij a timpului operator într-adevăr, ştim că valoarea medie a unei sume de variabile aleatoare este egală cu suma valorilor medii a fiecăreia dintre variabile : (2 39) (2 39) E{Xt + X2 + + Xn) = E(XJ + E(X2) + + E(Xn), (2 40) unde simbolul E reprezintă valoarea medie a variabilei aleatoare Xv Să reluăm exemplul din fig 2 22 presupunînd că ne găsim în cazul nesigur Valorile afectate arcelor din fig 2 22 vor reprezenta acum valorile cele mai probabile md Valorile atj şi bd sînt date în tabelul 2 8, unde am calculat şi prin formulele (2 37) şi (2 38) A fost reluat acelaşi exemplu din motive de simplificare a expunerii Trebuie să remarcăm totuşi că de fapt cele două probleme sînt diferite Vom reproduce pe fig 2 43 graful reprezen- tînd programul Fiecare arc reprezentînd o operaţie este afectat cu valoarea calculată şi care se poate citi pe tabelul 2 8 Numerele subliniate şi înscrise aproape de vîrfurile grafului reprezintă acum valorile medii ale datelor evenimentelor Et Astfel evenimentul E10 se va produce la o dată medie 47,99, evenimentul En la o dată medie de 41,16 şi evenimentul E12 care reprezintă terminarea lucrărilor la o dată medie de 60,82 ) Fig 2 43 Vom putea merge mai departe în cunoaşterea programului, ob- servînd că la capătul unui anumit număr de operaţii datele evenimentelor corespunzătoare sînt repartizate după legea normală, în virtutea teoremei limită centrale pe care o reamintim în cele ce urmează: Fiind dat un şir de variabile aleatoare independente Tabelul 2 8 p<> aH b<> mu 'u 5 11 8 8 1 Pi a 6 16 13 12,33 2,77 Pi e 6 12 9 9 1 3 6 4 4,16 0,25 Pt o 4 12 9 8,66 1,77 Pt e 3 10 6 6,16 1,36 ^3 4 5 12 7 7,50 1,36 ■^3 6 6 14 10 10 1,77 ■^s e 4 11 6 6,5 1,36 ■P3 3 7 14 9 9,5 1,36 ^4 7 7 16 10 10,5 2,25 ■P4 8 8 12 9 9,33 0,44 ^4 11 4 7 6 5,83 0,25 ^5 » 2 4 3 3 0,11 ■^e 7 4 7 5 5,16 0,25 ■P« s 2 4 3 3 0,11 ^e e 5 10 8 7,83 0,69 P7 11 2 5 4 3,83 0,25 ■^8,7 7 9 8 8 0,11 Pe » 3 6 4 4,16 0,25 Pe 10 2 6 5 4,66 0,44 P» n 7 16 13 12,50 2,25 PfhlO 3 6 5 4,83 0,25 Pion 8 17 13 12,83 2,25 Pll-10 5 9 6 6,33 0,44 Pun 9 19 17 16 2,77 cu dispersiile 2 2 2 tfl? • • • j dacă dispersiile sînt finite şi suma lor creşte indefinit cu n, dar astfel incit fiecare raport g?/S a? i tinde la 0, atunci variabila = Z„ - E(Zn) cu z, = (2,41) % urmează, la limită, o lege de repartiţie normală redusă, oricare ar fi distribuţiile variabilelor X{ Această înseamnă că o sumă de variabile aleatoare independente, fiecare dintre ele avînd o repartiţie oarecare, urmează, dacă numărul acestor variabile este suficient de mare, o lege normală (legea Gauss-Laplace), a cărei medie este suma mediilor variabile şi a cărei dispersie este suma dispersiilor variabile ) : E(Zn) = E(XJ + E(XJ + + E(XJ, (2 42) sau + • • • + (2-43) Teorema limită centrală se aplică bine, dacă graful programului este suficient de mare pentru evenimentele de la sfîrşitul programului şî în particular pentru evenimentul final, reprezentînd terminarea ansamblului lucrărilor Trebuie totuşi să remarcăm imediat că aceasta nu este valabilă decît cu condiţia ca timpii operatori să fie variabile aleatoare independente, ceea ce se admite în general, dar care poate să nu se întîmple în anumite cazuri In tot cazul, aceasta este o ipoteză care trebuie examinată în fiecare caz Astfel, în exemplul nostru, data evenimentului ZJ2 este repartizată după legea Gauss-Laplace de medie 60,82 şi de dispersie egală cu suma dispersiilor timpilor operatori care compun drumul critic Numind sau mai simplu Z12 data aleatoare a evenimentului E12, avem (^ ,) = °2 (*12) (ij a) + *2 +®2G4,b) + *2 (^ 11) + ^2( 7 la data 26,16 Această dată limită 26,16 ne convine, pentru că folosind-o nu riscăm să perturbăm prea mult data de realizare a ansamblului lucrărilor Bineînţeles, siguranţa ar fi mai mare luînd ca dată limită a lui Ee, de exemplu, valoarea 24 (fig 2 48) Sîntem mai siguri atunci că data lui E12 va fi dictată de drumul critic iniţial: (fig 2 43), mai curînd decît de operaţiile drumului (fig 2 45) Orice ar fi, ne vom conforma metodei PERT şi vom lua ca dată limită a evenimentului Es data 26,16 Astfel, în cazul general, data limită ti a unui eveniment se obţine scăzînd din data medie tn a evenimentului final En timpul mediu necesar realizării operaţiilor posterioare lui E{ Se procedează deci ca în 2 4, unde timpii operatori erau siguri, dar utilizînd de data aceasta timpii operatori medii (fig 2 43) : (2 44) dacă Tj n este timpul mediu necesar realizării operaţiilor posterioare lui E( în final, cunoaştem pentru fiecare eveniment data sa medie de realizare tt şi repartiţia corespunzătoare, la fel data sa limită Algoritmii folosiţi vor fi cei din 2 5, cu deosebirea că timpii operatori, care erau consideraţi determinaţi, sînt acum timpi operatori medii Vom defini un interval de fluctuaţie mediu (fig 2 49) (2-45) ti ] ti ]pentru evenimentul Ei şi o margine liberă medie pentru operaţia : h -ti- ta (2 46) Trebuie să remarcăm că se desemnează tot prin interval de fluctuaţie cantitatea ti — tj} cînd aceasta nu conduce la confuzii EEALIZAEEA UNUI CALENDAE Stabilirea unui calendar constă în fixarea datelor anumitor evenimente şi, în particular, a evenimentului final, corespunzător realizării programului Evenimentele interesante de luat în consideraţie sînt acelea ale drumului critic, pentru că ele condiţionează imperativ realizarea programului Cunoscînd distribuţia datei fiecărui eveniment de la sfîrşitul programului, este uşor de evaluat probabilitatea realizării sau a nerealizării acestor evenimente la data fixată în calendar, cu atît mai mult cu cît repartiţiile sînt, după cum s-a văzut înainte, normale sau cvasinormale Să reluăm exemplul din fig 2 43 şi să presupunem că este vorba de programul construcţiei unei vile şi că, din diverse motive, cum părătorul vilei a fixat arhitectului sau antreprenorului data 63 pentru livrarea cheilor în mînă Această clauză a fost adăugată la contract şi o depăşire a datei prevăzute generează pentru antreprenor o penalizare proporţională cu întîrzierea în livrare Acesta, înainte de a modifica cîte ceva la programul său, va evalua riscul de depăşire a scadenţei Dacă contractul prevede livrarea lucrărilor la data 63, se va evalua riscul unei întîrzieri în livrare : J?12 are o dată t12 repartizată după legea Gauss-Laplace cu valoarea medie 60,82 şi dispersia o?2 = 9,51 Variabila redusă 60,82 X == — ff12 urmează legea normală Pentru t12 = 63, x = urmează legea normală Pentru t12 = 63, x = redusă redusă = 0,706 « 0,7 3,085 = 0,706 « 0,7 3,085 şi se citeşte, pe tabelul funcţiei de repartiţie a legii normale reduse, probabilitatea corespunzătoare : 0,75 Astfel, probabilitatea realizării înainte de data 63 este 0,75; există deci 75 de ocazii favorabile Iar 100 ca vila să fie realizată înainte sau la data 63 Probabilitatea nerea- lizării la data 63 este 1—0,75 = 0,25 ; există deci 25 ocazii favorabile la 100 pentru ca vila să nu fie realizată la data 63 Riscul de a fi penalizat este deci departe de a fi neglijabil pentru antreprenorul care a acceptat această dată de livrare ) Dacă contractul prevede livrarea lucrărilor la data 64, probabilitatea nerealizării este 1 — 0,85 = 0,15 Dacă contractul prevede livrarea lucrărilor la data 65, probabilitatea nerealizării este 1—0,91 = 0,09 în sfîrşit, dacă data pusă în contract este 66, probabilitatea nerealizării la data prevăzută nu este mai mare decît 1 — 0,95 = = 0,05 Astfel riscul de a nu onora contracul nu este pentru antreprenor decît de 5% Aceste rezultate sînt foarte interesante pentru antreprenor: — Dacă el are posibilitatea de a face să se accelereze lucrările, va examina, potrivit datei puse în contract, dacă este cazul de a face aceasta sau nu Va ţine cont în această decizie de faptul că orice accelerare a lucrărilor antrenează în general o creştere a costului acestora — Dacă nu are posibilitatea de a face să se accelereze lucrările, îi va fi posibil atunci să vadă într-un mod precis dacă are interesul să accepte sau să refuze contractul ţinînd seama de data fixată şi de preţul fixat prin convenţie pentru penalizare Figura 2 50 dă aspectul legii de probabilitate şi a funcţiei de repartiţie a legii normale reduse UTILIZAREA UNUI CALCULATOR ELECTRONIC Exemplele date în paragrafele precedente erau exemple foarte simple de programe care nu comportau decît un număr foarte mic de operaţii în practică, programele de cercetare sau de dezvoltare, mai ales cînd analiza este foarte fină, comportă zeci şi uneori sute de operaţii Se impune în acest caz folosirea unui calculator electronic Două aspecte merită a fi subliniate : evaluarea iniţială a unui program şi controlul unui program în curs de realizare EVALUAREA INIŢIALĂ A UNUI PROGRAM Metoda expusă permite, după cum am văzut, pornind de la un proiect de program, să obţinem indicaţiile următoare: data de realizare medie tt şi data limită t* ale diverselor evenimente ; intervalele medii de fluctuaţie ale evenimentelor şi marginile libere ale operaţiilor necritice; probabilitatea realizării unui calendar prestabilit şi aceasta pentru evenimentele de la sfîrşitul programului Dacă proiectul programului este satisfăcător, adică dacă probabilitatea de realizare a unui calendar prestabilit este suficient de mare, acest proiect este adoptat şi programul trece în execuţie în caz contrar se poate pune problema remanierii lui şi pentru aceasta se va căuta prin simulare să se reducă timpii operatori ai operaţiilor critice, alocîndu-le de exemplu resurse suplimentare, poate în detrimentul operaţiilor cu margini libere importante în sfîrşit, cu sau fără schimbări, va putea fi stabilit un program coerent răspunzînd bine diverselor imperative CONTROLUL UNUI PROGRAM ÎN CURS DE REALIZARE Pe măsura avansării programului datele reale înlocuiesc estima- ţiile mai mult sau mai puţin precise făcute înainte Datele noi se introduc în calculator, ceea ce permite un control permanent al programului, precizînd evoluţia sa viitoare şi ţinînd cont de desfăşurarea sa reală Astfel se poate, pentru respectarea termenelor impuse de un contract, să se modifice sau să se accelereze anumite operaţii Responsabilul execuţiei programului are deci în mînă un instrument de previziune şi de control care îi permite să ia în timpul dorit, cînd este necesar, toate deciziile de corectare care se impun (fig 2 51) Fig 2 51 Forma documentelor de intrare şi de ieşire poate fi foarte variată şi poate să se adapteze, eventual, la necesităţile particulare ale fiecărui tip de program Luînd de o parte tabelul complet al estimărilor (tabelul 2 9), care va fi într-o oarecare măsură documentul de lucru de bază al metodei PERT, se poate examina producţia, automată a diverselor documente sau grafice : lista evenimentelor sau a operaţiilor critice, clasamentul evenimentelor sau al operaţiilor în funcţie de intervalul de fluctuaţie sau de marginea liberă; un grafic cu abaterile (în plus sau în minus) ale datelor evenimentelor calculate, în raport cu datele corespunzătoare fixate în registrul de scadenţă dorit sau calculat anterior; un grafic indicînd vizual progresul prevăzut al programului : datele de scadenţă calculate ale evenimentelor puse împreună cu intervalul de fluctuaţie corespunzător; graficul lui Gantt al operaţiilor programului indicînd în mod egal intervalele de fluctuaţie în sfîrşit, este important de subliniat că programul-maşină, element central al metodei, este făcut o dată pentru totdeauna, permiţînd adaptările eventuale pentru informaţiile de ieşire (natura» şi forma documentelor dorite, prezentarea etc ) în cazul în care programul studiat comportă un număr relativ restrîns de operaţii elementare, metoda poate fi aplicată cu rezultate bune fără ajutorul unui calculator, cel puţin pentru activitatea de evaluare iniţială a programului Pentru funcţia control al programului din contră, tratamentul manual al datelor devine repede nepracticabil dacă se urmăreşte efectuarea controalelor frecvente şi rapide Tabelul 2 9 Eveniment nr Data de realizare aşteptat! Data limită de realizare Interval de fluctuaţie mediu h Scadenta iniţială Probabilitatea scadentei Medie Dispersie Medie Dispersie 50 92 38 92 0 0 82 0,05 51 85 35 85 3 0 77 0,09 54 74 29 82 4 8 73 0,42 52 47 25 74 7 27 70 1,0 53 70 31 70 7 0 60 0,02 57 56 23 64 10 8 55 0,42 E1 (punerea în execuţie a programului) Tabelul 2 9, „Informaţii de ieşire”, corespunde unei părţi a grafului din fig 2 52 Fig 2 52 în concluzie, meritul esenţial al metodei PERT şi al variantelor sale constă în posibilitatea de a-ţi da seama, pe măsură ce se desfăşoară lucrările, de deplasarea urgenţelor şi de punctele unde trebuie făcut un efort deosebit De fapt, drumul critic nu este întotdeauna acelaşi El poate evolua în funcţie de fapte imprevizibile Rămîne întotdeauna posibil să se examineze repercusiunile lor asupra datei probabile de terminare a lucrărilor, să se introducă sarcini noi (sub- tratate, de exemplu) şi să se evalueze incidenţa lor în cîteva cuvinte, este vorba de un instrument îmbunătăţit şi perfecţionat permanent, care permite orientarea acţiunii la fiecare moment, gradînd eforturile GENERALIZAREA METODEI PERT INTRODUCERE în acest capitol se va prezenta o generalizare a metodei PERT Modelul este acum un model stochastic, aplicîndu-se bine in special programelor de cercetare, pentru care nu este întotdeauna posibil să se stabilească un graf bine definit EVENIMENTE DECIZIONALE ) în multe programe de cercetare sau de dezvoltare anumite evenimente au o semnificaţie deosebită, in sensul că de rezultatul lor depinde desfăşurarea ulterioară a programului Altfel spus, nu se poate stabili un program, pentru că de rezultatul anumitor operaţii sau al anumitor evenimente depinde conducerea proiectului Astfel de evenimente se vor numi decizionale Se pot găsi mai multe dintre acestea în acelaşi program Evenimentele decizionale intervin constant în programele de cercetare unde experienţele cerute depind de rezultatele anumitora dintre ele Ele se întîlnesc de asemenea în programele de dezvoltare Un exemplu simplu va permite cititorului de a situa mai bine problema Serviciul geologic al unei companii miniere, după o serie de măsurări geofizice pe teren, a localizat un zăcămînt important de plumb într-o regiune determinată din C^vennes Inginerul geolog estimează tonajul zăcămintului la 2 000 000 tone şi conţinutul său mediu, după eşantioanele găsite şi rezultatele cîtorva sondaje efectuate, la 2,5% Pb Direcţia Companiei este interesată de acest zăcămînt şi datorită optimismului geologului, decide crearea puţurilor de extracţie şi, eventual, a unei instalaţiii de spălare Se stabileşte atunci un program de lucrări în care este prevăzută, evident, o ■campanie de sondaje sistematice De rezultatele acestei campanii vor depinde operaţiile ulterioare ale programului într-adevăr, importanţa şi alegerea instalaţiilor miniere vor fi diferite, după cum tonajul şi calitatea zăcămîntului răspund sau nu speranţelor ivite Se va putea admite chiar că stabilirea unei instalaţii de spălare nu este necesară dacă importanţa zăcămîntului este prea mică; minereul va fi transportat atunci cu camioane şi pe calea ferată la uzina cea mai apropiată a companiei, care îşi va asuma sarcina tratării sale Din acest exemplu se vede că nu există numai un program ■de lucrări, ci mai multe, alegerea unuia dintre ele depinzînd de rezultatele campaniei de sondaje Rezultatul sondajelor este evenimentul •decizional, în acest exemplu Pentru simplificare, să presupunem că nu există decît trei eventualităţi referitoare la rezultatele campaniei de sondaje : zăcămînt foarte important justificînd construcţia a trei puţuri (dintre care două de aerisire şi de coborîre a personalului şi unul de ■extracţie) şi stabilirea unei instalaţii de spălare; zăcămînt de importanţă mijlocie care nu necesită decît două puţuri; zăcămînt puţin important necesitînd numai două puţuri şi mijloace de exploatare reduse Operaţiile necesare punerii în exploatare sînt următoarele : obţinerea permisului de exploatare ; stabilirea unei piste de 6 km; transportul şi instalarea la baza construcţiei a două sondeze; crearea clădirilor provizorii pentru biroul de planuri şi locuinţa muncitorilor sondori; gudronarea pistei; /) aducerea apei; g) campania de sondaje; h,) săparea şi echiparea a trei puţuri; h2) săparea şi echiparea a două puţuri; ij) transportul şi instalarea la bază a materialului de exploatare tip 1; transportul şi instalarea la bază a materialului de exploatare tip 2; transportul şi instalarea la bază a materialului de exploatare tip 3; h) construcţia birourilor şi a locuinţelor muncitorilor şi inginerilor ; j2) construcţia birourilor şi a locuinţelor muncitorilor şi inginerilor ; j3) construcţia birourilor şi a locuinţelor muncitorilor şi inginerilor ; k) trasări şi amenajarea fondului; construcţia unei instalaţii de spălare Obiectivul programului este punerea în exploatare a zăcămîn- tului Traseurile invocate în operaţia k sînt traseuri strict necesare începerii exploatării In sfîrşit, operaţia „Construcţia de locuinţe” pentru muncitori nu este de aceeaşi importanţă cu a reşedinţei miniere Vom reuni cele trei programe de lucrări diferite într-unul singur Obţinem atunci graful generalizat sau graful stochastie din fig 3 1 Operaţia c a fost fracţionată în c' şi c", deoarece se presupune că nu se poate întreprinde gudronajul drumului înaintea trecerii son- dezelor Operaţiile virtuale, punctate pe graf, au fost introduse res- pectînd regulile stabilite în capitolul 2 Evenimentul 7, care este eveniment decizional al acestui graf stochastie, este reprezentat printr-un romb Există pe acest graf trei evenimente terminale : Ew Eie şi Ea Este interesant acum să căutăm probabilităţile respective ale celor trei evenimente terminale Dacă nu se predă nici o informaţie suplimentară, se poate admite arbitrar că probabilităţile lor respective vor fi 0,33 Dar inginerul geolog poate avea destule motive pentru a crede că se găseşte foarte probabil în eventualitatea 1 I se va atribui atunci evenimentului En o probabilitate de 0,75, evenimentului ^18 probabilitatea 0,10 şi evenimentului Eu probabilitatea 0,15 Astfel, planificatorul va putea spune că există 75 cazuri favorabile din 100 de a se găsi în eventualitatea 1 şi de a avea o durată medie de lucrări corespunzînd evenimentului E21 Acest exemplu simplu, neavînd de altfel decît un raport îndepărtat cu problemele ridicate de punerea reală în exploatare a unui zăcămînt, a permis cititorului să înţeleagă în ce constă un eveniment decizional şi cum 6e poate el reprezenta pe un graf de program CÎTEVA PRECIZĂRI ASUPRA TERMINOLOGIEI FOLOSITE Introducerea noţiunii de eveniment decizional permite luarea în consideraţie a alternativelor care se prezintă în desfăşurarea anumitor faze ale unui program, multe dintre activităţi fiind cercetări de existenţă şi căutînd să demonstreze că anumite relaţii există sau nu Două drumuri care pleacă dintr-un eveniment decizional se vor numai conjunctive, dacă trebuie efectuate amîndouă operaţiile pe care ele le reprezintă Dacă trebuie realizată numai una sau alta dintrele ele, exclusiv, ele se vor numi disjunctive (fig 3 2) Arcele sau drumurile 1 şi 2 sînt aici disjunctive Pnjgmm PERT Pnjgmm PERT Fig 3 2 Graful de tip PERT introdus în capitolul 2 poate fi considerat ca un caz particular al grafului stochastic mai general, deoarece graful PERT nu conţine drumuri conjunctive Două configuraţii speciale ale grafului stochastic merită puţină atenţie Acestea sînt: — dependenţa între drumurile conjunctive; — confluenţele drumurilor conjunctive Dependenţa între drumuri conjunctive are loc cînd alegerea unei operaţii de efectuat depinde de răspunsul unui eveniment decizional situat pe un drum conjunctiv Acest caz este ilustrat de fig 3 3 Drumurile 1 şi 2 sînt drumuri conjunctive care conduc la evenimentele decizionale B şi respectiv C Să presupunem acum, că dacă răspunsul lui B este negativ, operaţiile ulterioare depind de răspunsul lui C Un mijloc de a reprezenta aceasta este repetarea pe graf a evenimentului decizional C, ceea ce s-a făcut (fig 3 3) în Cj Trebuie observat că (7 şi sînt identice şi că Cj trebuie să se situeze la aceeaşi dată cu C, sau după el A doua configuraţie specială a grafului stochastic este aceea a confluenţei de drumuri conjunctive într-un vîrf al grafului ajung două sau mai multe drumuri conjunctive, după ce fiecare dintre ele a trecut printr-un eveniment decizional Evenimentul D furnizează un astfel de exemplu (fig 3 3) în acest exemplu, drumurile Fig 3 3 4 şi 5 ajung la D şi pentru caD să fie atins este necesar ca răspunsurile lui B şi C să fie ambele afirmative ENUMERAREA EVENIMENTELOR TERMINALE POSIBILE Deoarece graful stochastic permite reprezentarea alternativelor, este clar că nu toate evenimentele terminale sînt realizabile Pot apărea anumite restricţii Dacă simbolul © reprezintă operaţia sumă disjunctivă şi simbolul n operaţia intersecţie, expresia următoare indică, sub formă logică, mulţimea evenimentelor terminale posibile ale programului Fig 3 4 Pe acest graf sint reprezentate numai evenimentele decizionale prin excluderea evenimentelor normale de cercetări din fig 3 4 Evenimentele terminale posibile sînt: {[(A © B) © X] n (Y © E)} ® {(G® H)@F} x n y = c ®j) (3 1), obţinem evenimentele (3-1) sau Pornind de la expresia posibile : (A şi (A şi (B şi y) E) T) (B şi E) (X şi Y) = (C sau D) (X şi E) G H F (3-2) terminale (33) Dependenţa dintre drumurile conjunctive reprezentate prin evenimentele decizionale 2 şi 2r restrînge evenimentele finale posibile, pe anumite dintre ele făcîndu-le imposibile Observînd că (X şi T) este egal cu (C sau D), lista definitivă a evenimentelor terminale posibile este atunci: A şi E B şi I C (3 4) (3 4)D X şi E G Acum, după ce toate evenimentele posibile au fost enumerate, se va căuta clasamentul acestora în termeni de probabilităţi de realizare PROBABILITATEA EVENIMENTELOR TERMINALE Probabilităţi apriori sînt acum afectate fiecăreia din alternativele evenimentelor decizionale Este ceea ce am făcut în exemplul din 3 2, atribuind eventualităţii 1 probabilitatea apriorică 0,75 şi eventualităţilor 2 şi 3 probabilităţile 0,10 şi respectiv 0,15 în legătură cu aceasta, trebuie remarcat că un eveniment decizional poate conduce la mai multe eventualităţi, trei în exemplul din 3 2, dar că este întotdeauna posibil să se introducă în graf numai evenimente decizionale cu două eventualităţi, ca în fig 3 4 Graful din fig 3 4 devine astfel acela din fig 3 5 : Eveniment final D X şi E F B şi Y H C A şi E G Eveniment final D X şi E F B şi Y H C A şi E G Probabilitatea 0,2058 0,1960 0,15 0,126 0,09 0,0882 0,084 0,06 Probabilitatea 0,2058 0,1960 0,15 0,126 0,09 0,0882 0,084 0,06 Entropia H = 2,88 Bţnox = 3 E = 0,96 Entropia H = 2,88 Bţnox = 3 E = 0,96 După cum s-a indicat mai înainte, graful de tip PERT apare ca un caz particular al grafului stochastic în care unui drum disjunctiv sau unei mulţimi de drumuri conjunctive, pentru fiecare eveniment decizional, îi este afectată o probabilitate egală cu 1 Probabilităţile evenimentelor terminale ale grafului din fig 3 5 pot fi calculate acum cu ajutorul probabilităţilor apriorice ale rezultatelor evenimentelor decizionale (tabelul 3 1) Tabelul 3 1 Evenimente finale Prob&bllitâţl compuse Probabilişti ala evenimentelor terminale A şi E p(A)-p(B/A) = (0,4) (0,3) (0,7) (1) = 0,084 b şi y p(B)-p(Y/B) = (0,6) (0,3) (0,7) (1) = 0,126 c p(C)-p(X)-p(Y) = (0,3) (0,7) (0,7) (0,6) = 0,0882 D p(D)p(X)p(Y)= (0,7) (0,7) (0,7) (0,6) = 0,2058 X şl E p(X) -p(EIX) = (0,7) (0,7) (0,4) = 0,196 G p(G) = (0,4) (0,5) 0,3) = 0,06 H p(H) = (0,6) (0,5) (0,3) = 0,09 F p(F) = (0,5) (0,3) = 0,15 FUNCŢIA ENTROPIE Clasamentul evenimentelor terminale posibile în funcţie de probabilităţile lor de realizare este un rezultat foarte util în exemplul dat, evenimentul D are cea mai mare probabilitate şi evenimentul (Z şi E) se plasează al doilea Evenimentul cel mai puţin probabil este evenimentul G O utilizare mai bună a acestor probabilităţi poate fi făcută folosindu-le comparativ, decît încercînd să judecăm un eveniment terminal pe baza probabilităţii sale absolute, deoarece probabilităţile evenimentelor terminale descresc, în general, în timp ce lungimea grafului şi numărul de evenimente decizionale cresc, în legătură cu aceasta, se poate folosi funcţia entropie, bine cunoscută în fizică sau în teoria informaţiei, care dă o măsură cantitativă a gradului în care anumite evenimente terminale sînt mai favorabile, din punct de vedere probabilistic, decît altele Această funcţie este definită prin relaţia (3-5) (3-5) unde H este entropia, pt — probabilitatea apriorică a evenimentului terminal de ordinul i şi N — numărul evenimentelor terminale De exemplu, într-o situaţie de alegere dichotomică cu probabilităţi apriorice cu 1/2, H este egală cu 1, dacă se ia logaritmul în baza 2 Dacă probabilităţile sînt 1 şi 0, H este egală cu 0 Funcţia entropie are o valoare maximală egală cu log2N, dacă toate probabilităţile apriorice sînt egale între ele Ca măsură a incertitudinii relative într-un sistem, expresia entropie relativă a fost definită astfel: (3-6) Acest concept poate fi utilizat ca o măsură a gradului în care anumite evenimente terminale sînt mai probabile decît altele în graful din fig 3 5 avem : = log2 8=3, H = — Pi l°g2Pi = 2,88, E = -^1 = o,96 3 ’ METODA PERT GENERALIZATĂ Prin introducerea evenimentelor decizionale şi prin afectarea de probabilităţi apriorice rezultatelor fiecăruia dintre aceste evenimente, a putut fi elaborată o generalizare a metodei PERT, utilă în special în programele de cercetări Probabilităţile apriorice afectate evenimentelor decizionale pot fi obiective sau subiective în acest ultim caz, ele traduc presupunerile sau intuiţiile specialiştilor consultaţi Să reluăm exemplul din fig 3 5, unde timpii operatori medii sînt daţi de tabelul 3 2 Figura 3 6 dă o reprezentare a grafului, ţinînd Tabelul 3 2 Operaţia Timpul operator mediu Operaţia Timpul operator mediu ^0>l 4 P 2 B 9 ■^1 3 4 Pt C 5 ^1 2 11 1 9 ^*3 2, 0 P1 S 12 6 ps r 3 p 2 ,B 13 P&>6 7 p3 t 8 P 6,G 9 ■^2-4 6 10 seama de scara timpului Pot fi obţinute informaţii cum sînt drumul critic, probabilităţile de realizare ale anumitor date ale unui calendar stabilit în prealabil (ponderate prin probabilităţile de a lua anumite drumuri) etc , urmărindu-se preocupările privind fiecare proiect particular Cunoaşterea probabilităţilor evenimentelor terminale împreună cu cele ale datelor medii de realizare procură indicaţii asupra desfăşurării viitoare a proiectului Astfel, evenimentul final D este cel mai probabil şi drumurile conducînd la el sînt 0-l -3-4 şi 0-1-2-4 Timpul său de realizare este de 30 săptămîni De asemenea evenimentul (X şi E) are ca probabilitate de realizare 0,196 şi obiectivul E va fi atins la o dată medie 24 Trebuie observat că, o dată cunoscut rezultatul evenimentului 2 şi dacă el este astfel incit să conducă la E, operaţia P3i conducînd la X poate fi întreruptă, deoarece pentru atingerea evenimentului 4 este necesar ca ambele operaţii X şi Y să fie efectuate Să observăm, în sfîrşit, că evenimentul 2, este obligat să se plaseze la aceeaşi dată cu evenimentul 2 sau după el Entropia relativă poate fi utilizată ca o măsură, sub forma unui număr unic, după gradul căruia anumite evenimente terminale sînt mai probabile decît altele O entropie relativă egală cu 1 va sugera că toate evenimentele terminale au aceeaşi probabilitate (sînt echi- probabile), iar o entropie relativă egală cu zero va indica degenerarea grafului într-un graf de tip PERT, unde un eveniment final este absolut sigur de a fi realizat şi toate celelalte nu pot avea loc în concluzie, informaţiile, obţinute prin metoda PERT generalizată, permite planificatorului sau organizatorului să compare alternativele unui program pe baza a trei caracteristici: şansele lor relative de realizare, operaţiile care conduc la realizare şi datele lor medii respective Această metodă se aplică în special programelor de cercetări, dar se aplică la fel de bine şi anumitor programe de dezvoltare OPTIMIZAREA FUNCŢIEI ECONOMICE A COSTURILOR INTRODUCERE în capitolul 2 au fost obţinute diferite informaţii foarte interesante cu ajutorul metodei PERT Stabilirea drumului critic sau a drumurilor critice a permis localizarea operaţiilor-cheie care comandă desfăşurarea unui program şi, în particular, data sa medie de realizare Calculul intervalelor de fluctuaţie şi al marginilor libere operatoare a pus în evidenţă elasticitatea evenimentelor sau a operaţiilor necritice şi latitudinea pe care o putem avea în punerea lor în execuţie Aceste aspecte diferite se regăsesc în faza „control al programului în curs de execuţie” Atenţia planificatorului poate fi atrasă de noi operaţii critice neprevăzute în programul iniţial El ştie atunci în ce direcţie trebuie să-şi îndrepte efortul, pentru a respecta calendarul impus în general, timpii operatori, fie că sînt determinaţi sau sînt aleatori, variază în funcţie de costul operaţiei Accelerarea unei sarcini se traduce în general printr-o creştere a costului său Vrem să luăm acum în consideraţie acest aspect şi obiectul acestui capitol este de a arăta cum este posibil, în anumite condiţii, de a optimiza un program UNELE CONSIDERAŢII GENERALE Fiecare timp operator tif poate varia între două limite du şi Z>„, impuse de consideraţii tehnice sau economice Timpul operator corespunde duratei normale a operaţiei : este durata corespunzătoare costului minimal Timpul operator du corespunde, din contră, duratei operaţiei, cînd aceasta este accelerată la maximum, iar preţul corespunzător este maximal Notînd cu pv costul operaţiei Pfj pentru o durată tijf putem admite că = f este o curbă de un tip analog celei desenate în fig 4 1 Costul creşte dacă tif descreşte pînă îndw, dincolo de care nu se mai poate efectua lucrul Pare destul de verosimil ca întîi costul să treacă printr-un minim foarte plat, apoi să crească pentru condiţii anormale de lucru corespunzătoare la mijloace insuficiente, în ceea ce priveşte mîna de lucru sau materialul Dacă se ia ca timp operator al fiecărei operaţii a programului durata normală Z>o, durata de execuţie totală a programului va fi durata normală : programul stabilit va fi un program normal Din contră, dacă se accelerează la maximum fiecare operaţie şi dacă pentru fiecare dintre ele se ia ca timp operator dtj, programul obţinut va fi programul numit accelerat al ansamblului de lucrări de realizat Trebuie să observăm că ne aflăm aici în situaţia în care timpii operatori au valori determinate şi nealeatoare; altfel spus, pentru un cost dat, durata unei operaţii este bine determină şi aceasta, pentru fiecare operaţie, se poate determina printr-o curbă de tipul celei din fig 4 1 Fig 4 1 Vom pune acum în evidenţă două aspecte interesante de optimizare a unui program, pornind de la informaţiile date de metoda drumului critic DIMINUAREA COSTULUI TOTAL AL UNUI PROGRAM în acest paragraf nu este vorba de căutarea optimului unui program, ci numai de a arăta cum permit informaţiile date de metoda PERT să i se diminueze costul total, acţionînd asupra operaţiilor necritice (4-1) (4-1)Să considerăm deci, pentru un ansamblu de lucrări, un program oarecare, unde timpii operatori au fost totuşi aleşi astfel încît Calculul drumului critic dă durata de execuţie a programului Costul total al programului astfel stabilit este suma costurilor tuturor operaţiilor: (4-2) unde piS = 0 dacă (i, j) e U (U reprezentând mulţimea arcelor care constituie graful programului) Pentru a diminua costul total al programului, durata de execuţie rămînînd aceeaşi, vor trebui micşorate marginile libere ale operaţiilor necritice în limitele impuse de restricţiile 4 1 într-adevăr, a diminua marginea liberă a unei operaţii înseamnă a mări timpul său operator, deci a diminua costul Făcînd aceasta s-ar putea întâmpla să se creeze noi drumuri critice Acesta va fi cazul cînd s-ar putea anula toate marginile libere ale operaţiilor care constituie un drum pornind de la la En Vom preciza toate acestea reluînd exemplul din 2 3 Timpii operatori tif sînt acum susceptibili să varieze între DiS şi dtj, mijlocind o creştere a costului fiecărei sarcini Valorile timpilor Dtj şi di} sînt date în tabelul 4 1 Vom presupune că creşterea costurilor este proporţională cu scăderea timpilor operatori Aceasta duce la presupunerea că se reduc curbele de tipul celei din fig 4 1, pentru fiecare operaţie, la segmente de dreaptă Coeficienţii de proporţionalitate c 23 = 6, durata lui P2|3 va putea fi lungită şi va fi de 5 săptămîni P215 avînd o margine liberă nulă, nu poate fi lungită P2,6 avînd o margine liberă 9 şi Z>2& fiind egală cu 7, durata sa va putea fi lungită şi va fi de 7 săptămîni etc Fig 4 2 Noii timpi operatori ai operaţiilor necritice vor fi t'l 2 = 8 «4 7 = 11 «8 10 = 6 «i 3 = 13 «4 8 = 9 «8 11 = 13 «l',4 = 10 «4 11 = 8 ^9 10 = 6 «2,3 = 5 «5 0 — 4 «10 12 = 13 «2,5 = 9 «« 7 = 7 «11 10 = 6 «2,6 = 7 «6 8 = 6 «11 12 = 18 «3 4 = 7 «8 9 = 8 «3 6 = 10 «7 11 = 5 «3 8 = 8 «8 7 = 8 «3 9 = 10 «8 9 = 4 Aceşti noi timpi operatori sînt sublimaţi Graful noului program va fi cel din fig 4 3 Anumite margini libere sînt anulate şi au apărut astfel noi drumuri critice Drumurile critice sînt: [1, 2, 3, 6, 8, 11, 10, 12 | [1, 3, 6, 8, 7, 11, 10, 121 , Noile margini din fig 4 3 sînt puse în paranteze Economia asupra fiecărei operaţii s-a obţinut făcîndu-se produsul lui co cu creşterea duratei operatorii Economia realizată asupra lui P1A este egală cu 200; asupra lui P2>3 este 300; asupra lui P2 6 este 150 ; asupra lu i P3 8 este 200 x 2 = 400; asupra lui P3i este 150; asupra lui P4>7 este 250; asupra lui PiM este 150 x 2 = 300 ; asupra lui Ps 9 este 400; asupra lui P6>7 este 150 x 2 = 300 ; asupra lui P8 8 este egală cu 300 x 3 = 900; asupra lui P7 u este 150 ; asupra lui P8-10 este 300; asupra lui P910 este 400 ; asupra lui Pu,12 este 500 Economia totală este suma economiilor realizate asupra fiecărei operaţii, adică 4700 F şi costul total al noului program din fig 4 3 este Nu trebuie însă să credem că s-a atins optimul pentru o durată de realizare de 61 săptămîni Consideraţiile precedente au căutat să arate pur şi simplu cum informaţiile date de metoda PERT permit ameliorarea unui program, prelungind operaţiile necritice atît cît este posibil ACCELERAREA UNUI PROGRAM LA CEL MAI MIC COST) Pentru a diminua timpul de execuţie total al unui program, trebuie diminuată durata unei operaţii critice Durata de execuţie totală va fi de asemenea redusă Dacă alegem operaţia critică care pentru o aceeaşi diminuare a timpului propune cea mai mică creştere a costului, vom accelera programul cu cele mai mici cheltuieli Şi aici, metoda folosită va fi intuitivă şi destul de banală Obiectul acestui paragraf este, ca şi al paragrafului 4 3, să arate importanţa economică a rezultatelor metodei PERT : drum critic şi margini libere Să presupunem, de exemplu, că avem de construit şi echipat un ansamblu hidroelectric Să luăm ca unitate de timp luna Operaţiile de realizat sînt: construcţia drumurilor de acces la uzină şi la carierele de exploatat pentru obţinerea materialelor (durata probabilă ta= 4); pregătirea carierelor de exploatare şi a fundaţiilor = 6); construcţia unui centru pentru administraţie şi personal, în general (te = 4); comanda materialului electric şi hidraulic (generatoare, turbine, conductă forţată etc ) (td = 12); construcţia uzinei (te = 10); /) construcţia barajului, a digurilor şi a deversorului de suprafaţă (tf = 24); construcţia galeriilor şi conductelor de aducţiune (is = 7); montajul uzinei şi al conductelor (tA = 10); t) încercarea funcţionării (t{ == 3) Graful acestui program este cel din fig 4 4 Evenimentele sau etapele sînt: pornirea lucrărilor; terminarea drumurilor de acces; terminafrea centrului şi a fundaţiilor; terminarea uzinei şi a conductelor de aducţiune ; terminarea ansamblului energetic; livrarea lucrărilor Fig 4 4 Am presupus că programul iniţial normal stabilit se ridică la 350 milioane F repartizaţi după cum se indică în tabelul 4 2 Pe acest tabel sînt indicate de asemenea duratele şi costurile corespunzătoare programului accelerat, unde toate operaţiile sînt reduse la minimul compatibil cu cerinţele tehnice Tabelul 4 2 Barclna Program normal Program aocelerat Costul accelerării pe lună Durata probabilă (luni) Cost probabil (milioane F) Durata probabilă Coet probabil a 4 5 2 15 5 b 6 11 5 30 19 c 4 3 2 11 4 d 12 150 9 180 10 e 10 10 8 20 5 f 24 147 19 212 13 9 7 18 6 30 12 h 10 4 7 25 7 i 3 2 350 2 5 528 3 Dacă vrem să reducem durata totală a execuţiei la 9 luni, trebuie să atacăm mai întîi sarcinile critice, adică sarcinile a, b, f şi i în mod precis, reducerea cu o lună a lui i va fi mai puţin scumpă (3 milioane), deoarece reducerea sarcinilor a, b şi f costă 5, 19 şi respectiv 13 milioane F Deci, vom urma această cale (fig 4 5) Nu este necesar să accelerăm sarcinile necritice; s-ar cheltui bani în Fig 4 5 pierdere; autorii americani o numesc eroare de alegere „To buy floating time” *) Să mergem mai departe şi să încercăm să cîştigăm o lună în plus, în acest caz, reducerea cu o lună a lui a ne va da costul cel mai mic, deoarece după cum se vede în tabelul 4 2, nu se poate continua cu reducerea lui i Se va continua în acelaşi mod accelerînd sarcina a, apoi sarcina / după cum se poate vedea în tabelul 4 3 Reducerea sarcinii f de la 21 la 20 luni va face să apară două noi sarcini critice, anume : e să fie P34 şi h să fie P4S (fig 4 6) în sfîrşit, dacă vrem să cîş tigăm încă o lună, adică să reducem durata totală de la 30 la 29 luni, putem realiza aceasta prin una din următoarele căi: (1 lună asupra lui f: cost 13 a cîştiga | asupra lui e : coBt 5 k Total: 18 a cîştiga { 1 lună asupra lui b : cost 19 Total: 19 Se vede că este avantajos să cîştigăm timp asupra lui / şi e Trebuie deci 18 milioane în plus faţă de situaţia de mai înainte pentru a cîştiga încă o lună în definitiv, alegerea va putea fi exercitată asupra programelor următoare (vezi tabelul 4 3) ') A cumpăra timp pierdut (N T ) Tabelul 4 3 Durata programului Costul (mllloa* ne F) Durata operaţiilor (luni) a b c | d ‘ I f 0 h 37 luni 350 4 6 4 12 10 24 7 10 3 36 - 353 4 6 4 12 10 24 7 10 2 35 - 358 3 6 4 12 10 24 7 10 2 34 - 363 2 6 4 12 10 24 7 10 2 33 - 376 2 6 4 12 10 23 7 10 2 32 - 389 2 6 4 12 10 22 7 10 2 31 - 402 2 6 4 12 10 21 7 10 2 30 - 415 2 6 4 12 10 20 7 10 2 29 - 433 2 6 4 12 9 19 7 10 2 28 - 452 2 5 4 12 9 19 7 10 2 VARIAŢIA ARBITRARĂ A COSTULUI OPERATOR ÎN FUNCŢIE DE( DURATĂ Este cazul foarte general cînd fiecărei operaţii i se poate asocia o curbă de cost în funcţie de durată, de tipul celei din fig 4 1 Fiecare timp operator variază între două limite : durata normală Dif corespunzătoare costului minim şi durata accelerată d{} Dacă într-un program se ia pentru fiecare operaţie durata sa normală, corespunzătoare costului minimal Dtj, se obţine un program normal al cărui cost total este el însuşi minimal Durata sa este Dn dacă n este indicele evenimentului final Dacă se ia pentru fiecare operaţie durata sa accelerată dijv se obţine un program accelerat Acesta nu este în general un program optimal, adică cu cost minimal, căci i se poate reduce costul prelungind operaţiile necritice, cum s-a indicat în 4 3 Din contră, durata sa de execuţie este minimală O vom numi dn dacă En este evenimentul final Pentru valorile situate între dn şi D„ există o infinitate de soluţii posibile Dacă reprezentăm variaţia costului total p în funcţie de durata de execuţie totală a ansamblului lucrărilor, obţinem o figură OABC care include toate soluţiile posibile, OA fiind locul soluţiilor optimale (fig 4 7) în cazul general, cînd curbele : costul în funcţie de durată sînt oarecare, metodele de cercetare operaţională, în stadiul cunoştinţelor actuale, nu par să fie capabile să determine curba OA, locul soluţiilor optimale Totuşi, metoda indicată în 4 3 permite determinarea punctului A Pornind de la un program unde toţi timpii operatori sînt putem reduce costul acestui program, micşorînd marginile libere ale operaţiilor necritice, în limitele permise Pe fig 4 7 putem traduce aceasta spunînd că evoluăm pe verticala ABC a lui Plecînd din punctul C, atingem A reducînd la maximum marginile libere ale operaţiilor necritice Plecînd de la programul corespunzător punctului 0 (timpi operatori Ba), putem diminua durata totală de execuţie reducînd arbitrar una sau mai multe operaţii critice şi continuînd astfel din aproape în aproape, observînd că, în general, se creează noi drumuri critice Evoluăm atunci pe figura OAB, care reprezintă mulţimea soluţiilor posibile, nescurtînd decît lucrările drumurilor critice Această abordare a problemei poate fi interesantă dacă figura OAB este suficient de aplatizată Poate fi construit un algoritm pentru a face metoda mecanizabilă Dacă acum, la fiecare etapă alegem operaţia sau operaţiile critice, care dau pentru o aceeaşi micşorare de timp cea mai mică creştere a costului, cum am făcut în 4 4, atunci evoluăm pe curba OA, locul soluţiilor optimale Din nefericire se înţelege că o astfel de metodă este inaplicabilă în cazul general cînd curbele cost-durată sînt oarecare şi cînd programul cuprinde mai multe sute de operaţii Ea nu se aplică uşor decît atunci cînd graful este foarte simplu (de ordinul a zece operaţii) sau atunci cînd, pentru fiecare operaţie, curbele pi} = au forme particulare (parabole, drepte, curbe convexe formate din segmente de dreaptă) Aceste cazuri particulare formează obiectul paragrafelor următoare VARIAŢIA LINIARĂ A COSTULUI ÎN FUNCŢIE DE DURATĂ PROGRAM LINIAR PARAMETRIC Se presupune că pentru fiecare operaţie Po curba costului în funcţie de durată, fie pij=f(tij), este o dreaptă (fig 4 8), adică: Pa = ~ ca hi (4 3) cu C fu C (4 4) Să considerăm un anumit program, a cărui durată de execuţie este tnf En fiind evenimentul final Ne propunem să minimizăm costul total al acestui program Aceasta s& traduce prin programul liniar / dij iii / X hj f»» pentru orice p e pB, j Meu (4 5) [MIN] £ (^ Ci /J, l (» )> eu unde ptt este mulţimea drumurilor mergînd de la vîrful 1 la vîrful n; U — mulţimea arcelor care constituie graful programului Fă- cînd să varieze durata totală de execuţie a programului tn şi înlocuind MIN prin MAX, obţinem programul liniar parametric, unde tn = X este parametrul care variază : , £ ttj X, pentru orice p e p „, [MAX] cij t,j (i i)ev (4-6) (4-6) dreaptă Metoda de rezolvare dreaptă Metoda de rezolvare ideile enunţate pe scurt în ideile enunţate pe scurt înîn cazul particular studiat, curba OA din fig 4 7 este o curbă convexă formată din segmente de clasică a programelor liniare parametrice ne permite să obţinem vîrfurile acestei curbe (fig 4 9) Cunoscînd curba OA, locul soluţiilor optimale, planificatorul poate alege programul care răspunde cel mai bine preocupărilor sale, îndreptîndu-şi atenţia fie asupra duratei totale de execuţie a ansamblului lucrărilor, fie asupra costului său total pe peSă observăm că metoda fo- iosi tă în rezolvarea unui program liniar parametric se bazează chiar 4 5 De fapt, pentru problemele care ne interesează şi anume optimizarea unui program, au fost elaboraţi mai mulţi algoritmi, care ţin seama de structura foarte particulară a programului liniar parametric (4 6) Aceştia sînt algoritmul lui Kelley şi algoritmul lui Fulkerson [231 Aceşti doi algoritmi sînt sensibil echivalenţi şi ambii au ca punct de plecare dualul programului (4 6) în paragrafele următoare vom expune complet algoritmul lui Fulkerson ALGORITMUL LUI FULKERSON ASPECTUL MATEMATIC AL ALGORITMULUI Notînd mereu cu t( data de realizare a evenimentului programul (4 6) se poate scrie sub forma : ti3 + h ~ 0 pentru orice (i, j) eU , (4 7) - «i+ tn eu Să transformăm variabilele /y în modul următor Se înlocuieşte fiecare f(j printr-o sumă de două variabile nenegative : fi) = fl fi i + fiji 2 > (4 19) noile variabile fiind supuse la restricţiile : fi i; i cu — S ^i} u,i) e u (i,D e c fi i; 1 — £ dij fi i; 2 a,i) e o Dacă punem Ci u t | Oii OO k = 1, k = 2; (4 22) A, du k = 1, k = 2 ; (4 23) i = n, i = 1, i = n, i = 1, (4 24) (4 24)programul dual se scrie neglijînd termenul constant: S X + XX- fa- k a i) ei'» /c ii,)) e v» o o (4-26) (4-26)(4 25) Acest program este acela al unei probleme de flux intr-o reţea dedusă din graful proiectului dedublînd toate arcele (fig 4 11) Fiecare arc (ij;&) al noului graf are o capacitate maximală Cu » şi îi este afectat un cost Sîntem astfel conduşi la o problemă de flux maximal, unde funcţia economică de considerat este expresia (4 26) Astfel de probleme fac de cîţiva ani obiectul a numeroase studii, şi sînt cunoscute pentru aceste cazuri particulare metode de calcul superioare metodelor generale ale programării liniare parametrice Timpii operatori optimali tu ai proiectului, corespunzători datelor tt produse în cursul procedeului, sînt daţi simplu de : = min [Aj, h - tj (4 27) ASPECTUL TEORETIC AL ALGORITMULUI Procesul de bază utilizat în algoritmul lui Fulkerson este un proces de marcare a anumitor vîrfuri ale grafului în general, acest proces este căutarea sistematică a unui lanţ, avînd anumite proprietăţi, de la vîrful 1 la n Cuvîntul lanţ prin opoziţie cu cuvîntul drum înseamnă că anumite arce pot fi traversate în sens invers orientării lor proprii Altfel spus, un lanţ este un şir de arce de o orientare oarecare ) Procesul de marcare se începe cu un flux ■ k şi cu datele tf care satisfac condiţiile : ti = o, di j k + tf — tj 0 => fij- k = Ci jlk- Pentru a simplifica notaţiile, să punem du; k = di t ,k + tj — tj (4 28) (4 29) (4 30) (4 28) (4 29) (4 30)Pentru X = tn aceste proprietăţi înseamnă că fluxul / minimizează (4 26) Aceasta se deduce din teorema asupra dualităţii unui program liniar sau se poate vedea direct, observînd că : k (i /) e fjn-k, apoi 6ă trecem la vîrful j Procedeul general este atunci: dacă vîrful p este marcat [i, k, +, eP ], se adaugă e„ la/^;* ; dacă vîrful p este marcat [i,k, e„] se scade e„ din fp^k', în ambele cazuri se trece la vîrful i Se repetă operaţia pînă cînd vîrful 1 este atins, adică e„ a fost adăugat la un arc (1, j ; k) Să ştergem acum toate marcajele şi să reluăm faza A Ambele eventualităţi de mai sus sînt reprezentate, respectiv, pe figurile 4 14 şi 4 15 Fig 4 15 Faza c | Schimbarea datelor evenimentelor Vîrful n n-a putut fi marcat Să definim mulţimile de arce Aj şi A2: Ai = { (i, 3 ; *) 1 i marcat, j nemarcat, di i t 0} (4 40) Să calculăm 3j = MIN [ - dU:t], (4-41) 32 = MIN [ d^j; *], (4 42) Să luăm 8 = min (8n S2) (4-43) Să schimbăm datele vârfurilor i scăzind 8 din toţi t, corespunzători unui vîrf i nemarcat Să ştergem acum toate marcajele şi să reluăm faza A Un arc al mulţimii Ai este prezentat pe fig 4 16 şi un arc al mulţimii A2 este reprezentat pe fig 4 17 în concluzie, algoritmul generează fluxuri succesive şi date tt care satisfac proprietăţile (4 28) — (4 30) Curba OA, locul soluţiilor optimale, este liniară între două valori succesive ale lui X = tn generate de algoritm Costul total al ansamblului lucrărilor este dat, pentru fiecare valoare a lui X calculată, de relaţia p (X) = K — [X® — 5] £ diii;k (4 44) * (» !) e c unde K este constanta : K = E [^-^ ,^1- (4-45) ( f (t(i) (fig 4 18 şi fig 4 19) Fig 4 18 Fig 4 18 Fig 4 19 Fig 4 19 4 9 EXEMPLU *) DE APLICARE A ALGORITMULUI LUI FULKERSON Un exemplu simplu va permite cititorului să se familiarizeze cu acest algoritm, al cărui aspect teoretic poate părea greoi Se consideră un program de lucrări al cărui graf este acela din fig 4 20 Tabelul 4 4 dă valorile lui , Diit şi citj Tabelul 4 4 % D 8 6 8 5,00 (1 3 12 36 4,66 21 27 6,00 25 50 1,00 «1 5 33 36 1,00 ^5 6 60 80 2,00 12 22 1,30 «7 8 50 90 2,00 ♦) Am reluat un exemplu dat tn , dar calculul este făcut aici cu ajutorul algoritmului lui Fulkerson Pornirea Se calculează datele evenimentelor luînd pe Z),,- ca valoare a timpului operator prin metoda generală (prezentată în capitolul 1) şi ale cărei formule sînt amintite în (4 32) Se ia pentru fiecare arc un flux nul Programul este acela care corespunde punctului zero din fig 4 7 în fig 4 21 datele evenimentelor sînt încadrate Numerele din paranteză sînt marginile libere ale operaţiilor Se vede că Xj = = tn = 130 Procesul de afişaj 1° Pentru a întreprinde procesul de afişaj trebuie să începem prin a calcula pentru fiecare arc valorile di /; 1 = Ditj + ti — tt djj: 2 = di ) + ti — tj Aceste valori sînt reproduse în tabelul 4 5 Faza A în coloana di j;2 nu există zerouri; se trece atunci direct la partea a doua a procesului de afişaj, utilizînd coloana lui di , 1: di,2 i = 0 şi /i 2, i = 0 + li~ 5 = 1,00 ; vîrful 5 va fi deci marcat La fel vîrful 6 va fi marcat Begula (4 36) permite marcarea vîrfurilor 3 şi 2 etc Se obţine graful din fig 4 26 Vîrful 8 a fost marcat, se poate deci mări fluxul trecînd la faza B Faza B Fluxul este mărit cu o unitate în arcele: (6, 8; 1), (5, 6 ; 1) şi (1, 5 ; 1), el rămănînd neschimbat în celelalte arce Faza A Se reîncepe procesul de afişaj cu acest nou flux Nu se mai poate marca vîrful 5 pentru că, deşi = 0, avem acum = C1 5 = !• Singurul vîrf care se poate marca este vîrful 7 (fig 4 27) Vîrful 3 nefiind marcat se trece la faza C Faza C Se consideră numai mulţimea Ar di,2:2 = 3 dl,3; 2 = 52, — dl,3; 1 = 28, dl,6; 2 = 3 dî e; 2 = 52, — dî B: i = 12, 8, = MIN[ - dW: * ] = 3, 8 = 3 De aici rezultă noile date ale evenimentelor nemarcate din figura precedentă Noile date sînt încadrate înfig 4 28 Fluxurile nu figurează în acest graf Aplicarea formulei (4 44) dă p (121) = K - ■^121 unde Z124 este cantitatea 20 + 50 + 52 + 8 pe care am întîlnit-o în calculul lui Ap pentru X = 124 Creşterea costului total Ap în raport cu p (130) este Ap = p (121)-p (130) = 130 + 124 - 242 = 12 Mai general, se poate arăta că la fiecare iteraţie creşterea costului total Ap este dată de formula A p = Zo + (v — v0) Xo — v X, unde Zo, v0 şi Xo se referă la iteraţia precedentă şi Z = Zo + + (® — ®0) Xo La fel, creşterea costului în raport cu iteraţia precedentă este A'p = Ap — Ap0, în cazul de faţă, avem A'p = 12 — 6 = 6 Graful programului optimal ta = X = 121 este cel din fig 4 29, Numerele încadrate reprezintă datele evenimentelor Numerele înscrise pe arce sînt timpii operatori Numerele dintre paranteze sînt marginile libere ale operaţiilor corespunzătoare : P1>3 are ca margino liberă 25 ca şi P4>6; P7i8 are ca margine liberă 9 Există două drumuri critice : şi 3° înainte de a ne reîntoarce la faza A, să stabilim un nou tabel (tabelul 4 7), analog tabelului 4 5, cu datele calculate mai înainte Tabelul 4 7 Operaţia li — tj di'j; I = Di ) +(i - ai-j; 2 " i 1 + - 0 di ;a ■= dt j + fţ — h (1,2) 0 - 11 = - 11 20 - 11 = 9 11 - 11 = 0 (2,3) 11 - 52 = - 41 50 - 41 = 9 30 - 41 = - 11 (3,6) 52 - 104 = - 52 52 - 52 = 0 32 - 52 = - 20 (6,8) 104 - 112 = - 8 8 - 8 = 0 6 - 8 = - 2 (1,3) 0 - 52 = - 52 36 - 52 = - 16 12 - 52 = - 40 (2,4) 11 - 38 = - 27 27 - 27 = 0 21 - 27 = - 6 (4,6) 38 - 104 = - 66 50 — 66 = — 16 25 - 66 = - 41 (1,5) 0 - 33 = - 33 36 - 33 = 3 33 - 33 = 0 (5,6) 33 - 104 = - 71 80 - 71 = 9 60 - 71 = - 11 (1,7) 0 - 22 = - 22 22 - 22 = 0 12 - 22 = - 10 (7,8) 22 - 112 = - 90 90 - 90 = 0 50 - 90 = - 40 Faza A Există două zerouri în coloana di ,;2 • Vîrfurile 5 şi 2 vor fi deci marcate [1, 2, + , oo) dl,7; 1 = 0 şi /1,7; 1 = 0 7;1 = 0, avem acum /1;7,1 — Ci; 7 1 = Cj; 7 = 1,3 Vîrfurile marcate sînt indicate pe fig 4 36, ca şi noile fluxuri Vîrful 8 nefiind marcat, se trece la faza C Faza C Se consideră numai mulţimea Aj pentru că nu există arce (i, j -,k) cu i nemarcat şi j marcat dl 3; 2 = 40, — 21 3; i = 16, dl,7; 2 = 10, d2,3; 2 = 11, ds 6; 2 = 11, d4,#;2 1 41, d4,6;l = 16, = MIN [- 2,- ,;i] = 10, 8 = 10 De aici rezultă noile date ale evenimentelor nemarcate din figura precedentă Noile date sînt încadrate în fig 4 37 Fluxurile nu figurează în acest graf Creşterea costului total în raport cu p (130) este Ap = Z112 + (4,8 - 3,5) 112 - 4,8 -102 ■^102 sau Ap = 435,5 + 145,6 — 489,5 = 91,5 Creşterea costului în raport cu p (112) este A' p = 91,5 - 43,5 = 48 Timpii operatori sînt calculaţi prin formula (4 27) : = min [Dm, (t, — tJJ Graful programului optimal ts= X = 102 este cel din fig 4 38 Numerele încadrate reprezintă datele evenimentelor Numerele înscrise pe arce sînt timpii operatori, iar numerele din paranteze sînt marginile libere ale operaţiilor corespunzătoare -f*13 are ca margine liberă 6, la fel ca şi Pi 6 Există trei drumuri critice : ; şi 5° înainte de a reîncepe faza A, trebuie să stabilim un nou tabel (tabelul 4 9), analog tabelului (4 5), cu datele calculate mai înainte Tabelul 4 9 Operaţia ti ~ ti di fci =• Dij + ~~ tj dj ja “ d# + ti — bj 1,2 0 - 11 = - 11 20 - 11 = 9 11 - 11 = 0 2,3 11- 42 = - 31 50 - 31 = 19 30 - 31 = - 1 3,6 42 - 94 = - 52 52 - 52 = 0 32 - 52 = - 20 6,8 94 - 102 = - 8 8-8=0 6 - 8 = - 2 1,3 0 - 42 = - 42 36 - 42 = - 6 12 - 42 = - 30 2,4 11 - 38 = - 27 27 - 27 = 0 21 - 27 = - 6 4,6 38 - 94 = - 56 50 - 56 = - 6 25 - 56 = - 31 1,5 0 - 33 = - 33 36 - 33 - 3 33 - 33 = 0 5,6 33 - 94 = - 61 80 - 61 = 19 60 - 61 = - 1 1,7 0 - 12 = - 12 22 - 12 = 10 11 - 12 = 0 7,8 12 - 102 = - 90 90 - 90 = 0 50 - 90 = - 40 Faza A Există acum trei zerouri în coloana d!j ;2 Vîrfurile 5, 2 şi 7 vor fi deci marcate cu dî s: i = 0 şi /t 8: i = 1,30 (,, (t,- /,)] Graful programului optimal ta ■■= X = 99 este cel din fig 4 46 Pl 3 are ca margine liberă 5, la fel ca şi P4t9 Există trei drumuri critice : , şi 7° înainte de a trece din nou la faza A, trebuie stabilit un nou tabel (tabelul 4 11), analog tabelului 4 5, cu datele calculate mai înainte Tabelul 4 11 Oper&fl* ; 2 " ~ G' 1,2 0 - 11 = - 11 20 - 11 = 9 11 - 11 = 0 2,3 11 41 = - 30 50 - 30 = 20 30 - 30 = 0 3,6 41 - 93 = - 52 52 - 52 = 0 32 - 52 = - 22 6,8 93 - 99 = - 6 8-6=2 6 - 6 = 0 1,3 0 - 41 = - 41 36 - 41 = - 5 12 - 41 = - 29 2,4 11 - 38 = - 27 27 - 27 = 0 21 - 27 = - 6 4,6 38 - 93 = - 55 50 - 55 = - 5 25 - 55 = - 30 1,5 0 - 33 = - 33 36 - 33 = 3 33 - 33 = 0 5,6 33 - 93 = - 60 80 - 60 = 20 60 - 60 = 0 1,7 0 - 12 = - 12 22 - 12 = 10 12 - 12 = 0 7,8 12 - 99 = - 87 90 - 87 = 3 50 - 87 = - 37 Faza A Se atinge vîrful 8 prin arce (i,j :2) : drumul (fig 4 47) Calculele s-au terminat S-a atins punctul A pe curba soluţiilor optimale din fig 4 7 Tabelul 4 12 rezumă calculul paginilor precedente reprezintă marginea liberă a operaţiei Pi} sau (i, j)] Tabelul 4 12 Operaţiile Pii Timpii operatori: tţ j 0 1 2 3 6 (1,2) 20 14 11 11 11 11 11 (2,3) 50 50 50 41 31 30 30 (3,6) 52 52 52 52 52 52 52 (6,8) 8 8 8 8 8 8 6 (1,3) 36 36 36 36 36 36 36 ■^ij> 34 28 25 16 6 5 5 (2,4) 27 27 27 27 27 27 27 (4,6) 50 50 50 50 50 50 50 25 25 25 16 6 5 5 (1,5) 36 36 33 33 33 33 33 (5,6) 80 80 80 71 61 60 60 1^4 » 6 — — — — — — (1,7) 22 22 22 22 12 12 12 (7,8) 90 90 90 90 90 89 87 18 12 9 - — — — fg = X 130 124 121 112 102 101 99 A/q — 6 3 9 10 1 2 Ap — 6 12 43,5 91,5 97 111 A'p — 6 6 31,5 48 5,5 14 A'p/Af8 — 1 2 3,5 4,8 5,5 7,0 Se poate trasa de asemenea curba OA, locul soluţiilor optimale (fig 4 48) Aceasta este o curbă convexă formată din segmente de dreaptă Pantele succesive ale acestor segmente de dreaptă sînt date în ultimul rînd al tabelului 4 12 Acest exemplu de apbcare a algo- 12J + 8) Va*,*,*») Fig 4 48 ritmului lui Fulkerson a putut părea plictisitor Totuşi nu este inutil să tratăm complet un exemplu şi să studiem în amănunţime toate ■operaţiile pentru a avea siguranţa că algoritmul s-a înţeles bine Evident, îndată ce graful posedă mai mult de 15 arce, calculul devine imposibil de realizat cu mijloacele obişnuite şi trebuie încredinţat unui calculator electronic 4 10 PROGRAM OPTIMAL Cu ajutorul programării liniare parametrice (vezi 4 6) sau mai bine cu ajutorul algoritmului lui Fulkerson (vezi 4 7 şi 4 8) a fost posibil să se traseze curba O A, locul soluţiilor optimale Dar oprirea unei instalaţii sau întârzierea în livrarea unui ansamblu (uzină, vapor, imobil etc ) se poate traduce, în anumite cazuri, printr-un cost care este funcţie de durata totală a lucrărilor Acesta va fi, de exemplu, „lipsa de cîştig”, în cazul unei instalaţii industriale, în exemplul simplu al construcţiei unei vile, acest cost va fi pierderea pentru cumpărător cauzată de cheltuielile de hotel Daeă putem estima sau calcula acest cost, c = rw 0 (4 49) Dispersia este 0 şi B > 0, şi a cărei densitate de probabilitate este f (/) = 0, — oo unde ta este valoarea aleatoare a arcului at, adică, mai precis, durata aleatoare a operaţiei reprezentate prin arcul p(it7) (II 5) este repartiţia vectorului timp operator tv~ Pentru anumite valori / = ( este valoare medie a drumului critic mergînd de la vîrful 1 la vîrful i Dar valoarea lt (t) nu depinde de valorile arcelor fasciculelor pentru că nici un drum de la 1 la i nu utilizează vreunul din arcele p (/,) = P( max [ £ p (tv) + tf), (tv) (f{ + max (/, + tn ,fi + tf) (11 23) Inegalitatea (11 18) fiind adevărată pentru avem (11 24) (11 24)yi+1 > max (& + «i, • • - , ? max IE P (M [*i + AL • • • S p(to)max[ei + , et + tj (11 31) 'c Inegalitatea (11 25) fiind adevărată pînă la i şi folosind (11 16), avem ®i+i > £ P (M max LA + tlt ,fi + tj = /i+1 (11 32) •u Această relaţie verifică (11 25) şi demonstrează (11 17) Trebuie să observăm că inegalităţile gt lt, sînt valabile fie că fasciculele sînt independente sau nu Dimpotrivă, ipoteza de independenţă între fascicule este esenţială pentru (11 17) CÎTEVA EXEMPLE NUMERICE în acest paragraf prezentăm cîteva exemple numerice pentru a compara pe gn, f„ şi en în fiecare dintre aceste exemple am presupus aceeaşi repartiţie uniformă pentru fiecare arc şi independenţa între arce Numerele situate pe arce sînt valorile pe care le poate lua timpul operator corespunzător Pentru exemplele 2 şi 3, calculele au fost efectuate de un calculator electronic Exemplul 1 (fig II 10) Exemplul 2 (fig II 11) Exemplul 2 (fig II 11)= 3,00, /4 = 3,22, e4 = 3 32 >4 0 1 2 3 4 5 6 Frecvenţa 1 11 50 74 71 27 9 >4 0 1 2 3 4 5 6 Frecvenţa 1 11 50 74 71 27 9 Distribuţia valorilor drumului critic 3io — 2,50, f10 — 3,50, e10 — 3,95 Exemplul 3 (fig 11 12) 3e — 2,50, f6 — 2,92, ee — 2,95 Din aceste trei exemple se constată că aproximaţia f este mai bună decît aproximaţia g 111 SUBDIVIZIUNEA OPERAŢIILOR ŞI REALIZAREA UNUI CALENDAR*) în această anexă vom înfăţişa efectele subdiviziunii operaţiilor asupra probabilităţii de realizare a unui calendar prestabilit Mai prcis, problema tratată în această anexă este următoarea : Fiind dată o activitate sau un şir de activităţi care se termină cu un eveniment cu o dată fixată, să se determine efectul subdiviziunii unei activităţi (sau a mai multora) asupra probabilităţii calculate de realizare a evenimentului la data indicată Se presupune că o operaţie, pentru care s-au stabilit estimaţiile optimistă, cea mai probabilă şi pesimistă ale timpului operator, a, m şi b, respectiv, este divizată în n subactivităţi, conectate în serie, cu estimaţiile optimistă, cea mai probabilă şi pesimistă a;, wij şi b, ale timpului operator, unde i = 1,2, ,,n Se mai presupune că £ a, = a, £ bk = b, (III l) i = l i-1 i = l şi că timpii operatori ai subactivităţilor sînt variabile aleatoare independente : d = i = l unde d = b — a şi d, = bt — at, i = 1, n Utilizînd ecuaţiile repartiţiei |3, referitoare la valoarea medie şi la dispersia timpului operator, avem pentru activitatea originală, ♦) Aceasta expunere este extrasă din articolul înaintea subdiviziunii : E(t) = — (a + 4m+b), 6 E(t) = — (a + 4m+b), 6V (t) = [- (b - a)]2 = d2/3C> 6 După subdiviziune, avem : Ea (t) = £ 2? (t ) = A (flj + 4 + bt) t = l & i=L = — (a + 4m + b), 6 = — (a + 4m + b), 6F„W - t v W = £ [A (6 - a,) ]2 = dt Vn(t) (III 3a) (III 3b) (III 4a) (III 4b) (III 5) (III 4b) (III 6a) (III 3a) (III 3b) (III 4a) (III 4b) (III 5) (III 4b) (III 6a)Neconsiderînd decît cazurile de nedegenerare, scriem în metoda PERT probabilitatea de realizare a unui eveniment înainte de (sau la) data fixată în calendarul T„ se obţine calculînd TS — E(T) v = » KV(!C) (III 7a) unde E(T) şi V(T) sînt respectiv valoarea medie şi dispersia datei evenimentului luate din tabela legii normale (Gauss-Laplace) reduse, în această anexă, t reprezintă timpul operator al unei operaţii şi T este data de realizare efectivă a unui eveniment, T, fiind data prevăzută în calendar, e-4-d e-4-d (III 7b) (III 7b) Să presupunem acum că activitatea care a fost divizată se găseşte pe drumul folosit pentru calcularea lui E(T) şi V(T) Putem presupune de altfel că activitatea sudivizată este prima Dacă V(t) este dispersia timpului operator al primei operaţii, înaintea subdiviziunii, şi u2 este suma dispersiilor corespunzătoare operaţiilor următoare pe acest drum, atunci dispersia datei evenimentului este Să presupunem acum că activitatea care a fost divizată se găseşte pe drumul folosit pentru calcularea lui E(T) şi V(T) Putem presupune de altfel că activitatea sudivizată este prima Dacă V(t) este dispersia timpului operator al primei operaţii, înaintea subdiviziunii, şi u2 este suma dispersiilor corespunzătoare operaţiilor următoare pe acest drum, atunci dispersia datei evenimentului este V (T) = V(t) + V (T) = V(t) + (III 8) (III 8) Din ecuaţiile (III, 7 a) şi (III 8) se obţine Din ecuaţiile (III, 7 a) şi (III 8) se obţine Ts ~E(t) V V(0 + Ts ~E(t) V V(0 + (III 9a) (III 9a) înainte de subdiviziunea primei activităţi şi înainte de subdiviziunea primei activităţi şi Ts - E(t) y” ~ V Vn(0 + Ts - E(t) y” ~ V Vn(0 + (III 9b) (III 9b) după subdiviziune după subdiviziune Din ecuaţiile (III 6b) şi (III 9) se deduce: ly„l>lyl (in io) şi sînt posibile trei cazuri: pentru Ts > E (T), probabilitatea de a îndeplini evenimentul înainte de (sau la) data calendarului va creşte; pentru Ts :25 215 r:25 135 BIBLIOGRAFIE Este evident că din lipsă de spaţiu nu putem cita toate articolele publicate despre sistemul PERT Dăm aici numai articolele de bază pe care le-am utilizat la redactarea acestei broşuri şi pe care cititorul le va putea consulta cu interes Cărţi Behge, C Teoria grafurilor şi aplicaţiile ei (trad din 1 franceză) Bucureşti, Editura tehnică, 1969 |2] Kaufmann, A , Cruon, R La programmation dynamique Dunod, 1965 Kaufmann, A Metode şi modele ale cercetării operaţionale (trad din 1 franceză) Bucureşti, Editura ştiinţifică, 1967 Aceste lucrări nu se referă la metoda drumului critic, dar in ele se vor găsi diferite elemente ale teoriei grafurilor şi ale programării dinamice care slnt utilizate in mai multe paragrafe din această carte Articole Capitolul 2 P E R T , Pliase I & II Summary Report, Special Projects Office Dept of Navy, Washington D C Malcolm, D G , Roseboom, J H , Clark, C E , Fazar, W Application of a Technique for Research and Development Program Evalualion, Operations Re- search, 7, 5, 1959 Fuchs, C Une măthode d’Analyse et d’ăvaluation des programmes : le P E R T , Ghefs d’Entreprise, 1962, A Ftude Centre Franţais de R O , P E R T , 12/61 Aceste trei articole sint articolele de bază pentru capitolul 2 Se va putea consulta oricare din cele trei Grubbs, Frank E Attempts to validate certain P E R T statistics or «picking on P E R T • Opns Res , 10, 6, 1962 Clare, Charles E , « The P E R T model for distribution of an activity time» Opns Res , 10, 3, 1962 Articolele şi dau precizări asupra repartiţiei |J utilizate în sistemul PERT Fulkerson D R , Expected criticai path lengths in P E R T networks, Opns Res , 10, 6, 1962 Healy, Thomas L Activity Subdivision and P E R T probability statements Opns Res , 9, 3, 1961 Ventuba, E Impressions sur Ies tendances actuelles de la Recherche opdrationnelle aux U S A , Revue Franţaise de Recherche opărationnelle (SOFRO), 25 Roy, B Graphes et Ordonnancements SOFRO, 23 Articolele şi reprezintă o bună introducere în metoda PERT Roy, B Contribution de la thiorie des graphes ă l’ătude des problămes d’ordon- nancement 2e Confirence Internaţionale sur la R O , Aix-en-Provence, 1960 Roy, B , Contribution de la thiorie des graphes ă l’6tude de certains problfcmes- liniaires C R Acad Sciences, T 248, p 2 437, 1959 ,15] Roy, B , Cheminement et connexiU dans Ies graphes Application aux problfemes- d’ordonnancement METRA, sdrie speciale, n° 1, 1962 Roy, B , Simonnabd, M «Nouvelle măthode permettant d’explorer un ensemble des possibilitis et de diterminer un optimum » Revue Franţaise de Recherche opârationnelle, n 18 Algan, M , Roy, B , Simonnabd, M Principe d’une mithode d’exploration de certains domaines et application ă l’ordonnancement de la construction des grands ensembles Cahiers du Centre de Mathdmatiques et de Statistique appliquies aux Sciences sociales Ed Institut Solvay, Bruxelles, 3, (1962), p 41 Articolele — vor permite cititorului să aprofundeze aşa-zisa metodă a potenţialelor de care am vorbit la sfirşitul capitolului 2 Capitolul 3 Eisner, Howard A generalized network approach to the planning and scheduling of a research project Opns Res 10, 1 Information Scientifique Bull, La programmation dynamique, par A Kaufmann Capitolul 4 Cabb£ M^thode du Chemin Critique Service Technologie Shell Petit Couronner Rouen Kelley, James E Criticai Path planning and Scheduling — Mathematical Basis, Opns, Rs , 9, 3, 1961 Kelley, James E Jr , & Morgan, R Walker Criticai Path Planning and Scheduling, 1959 Proceeding of the Eastern Joint computer conference, p 160 — 173 Fulkerson D R A Network Flow computation for project cost curves — Rând Paper, P — 1947, 1960 — Management Science, 1961 Bildson, R A , Gillepsie, J R „Criticai Path planning P E R T integration” Opns Res 10, 6, 1962 Cititorul care doreşte o bibliografie foarte completă a articolelor americane apărute- despre metoda PERT va putea să consulte : Bigelow, Clifford George „Bibliography on Projet Planning and Control by Network Analysis” 1959 — 1961 Opns Res 10, 5 1962 DICŢIONAR DE TERMENI FOLOSIŢI ÎN LUCRARE Metoda drumului critic — (Măthode du chemin critique — Criticai Path Planning — Criticai Path Method) PERT — estimare şi revizuire a planurilor (Program Evaluation Research Task Program Evaluation and Review Technique) Opera[ii — lucrări, activităţi (Opărations, jobs, tasks, activities): Pit Eveniment — obiectiv (ăvănement, event): Program — (programme, program, planning) Calendar — (calendrier, schedule, scheduling) Craf — diagramă cu săgeţi (graphe, dirccted network, flow chart, linear graph) Timp operator — durata unei operaţii (temps opăratoire, job tome, duration of an acti- vity, elapsed time): tit Timp operator mediu — durata medie a unei operaţii (temps opăratoire moyen, expected elapsed time, mean time) : t^ Interval de fluctuafie — (intervalle de flottement, (floating time, slack) Data aşteptată a unui eveniment — data cea mai apropiată (date attendue d’un ăvănement; earliest time) : 1, Data limită a unui eveniment — data cea mai tlrzie (date limite d’un ăvănement, latcst time) : Arc — săgeată (arc, arrow) Vtrf — nod (sommet, node) Drum — (chemin path) Durata minimală — estimaţie optimistă (durăe minimale, optimistic time): av Durata maximală — estimaţie pesimistă (durăe maximale, pesimistic time) : b{, Durata cea mai probabilă — (dur^e la plus probable, more probable time): m^ CUPRINS Cuvint înainte Prefaţă GRAFURI DE ORDONANŢARE Introducere Noţiunea de grai Relaţia de ordine strictă intr-un graf conex fără circuite Descompunerea unui graf conex fără circuite in nivele Găsirea unei ordonanţări STABILIREA UNUI PROGRAM DE CERCETARE SAU DE PRODUCŢIE Introducere Reprezentarea unui program printr-un graf Data realizării ansamblului de lucrări Drumul critic Intervale de fluctuaţie Margini ale operaţiilor Căutarea drumului critic Algoritmii utilizaţi Două exemple de aplicaţie Operaţii cu durată aleatoare Realizarea unui calendar Utilizarea unui calculator electronic GENERALIZAREA METODEI PERT Introducere Evenimente decizionale Cîteva precizări asupra terminologiei folosite Enumerarea evenimentelor terminale posibile Probabilitatea evenimentelor terminale Funcţia entropie Metoda I’ERT generalizată OPTIMIZAREA FUNCŢIEI ECONOMICE A COSTURILOR Introducere Unele consideraţii generale Diminuarea costului total al unui program 108 Accelerarea unui program la cel mai mic cost 112 Variaţia arbitrară a costului operator în funcţie de durată 116 Variaţia liniară a costului in funcţie de durată Program liniar parametric 118 Algoritmul lui Fulkerson 120 Detalii ale procedeului iterativ utilizat in algoritmul lui Fulkerson 121 Exemplu de aplicare a algoritmului lui Fulkerson 129 Program optimal 154 Optimizarea unui program cind duratele operaţiilor sînt aleatoare 155 Concluzii 161 A X E X E I Repartiţia P 166 Valoarea medie a unui drum critic intr-un graf de program 168 Ipoteze şi notaţii 168 Aproximaţia ft a valorii unui drum critic 171 Citeva exemple numerice 175 Subsiviziunca operaţiilor şi realizarea unui calendar 177 Aplicarea metodei potenţialelor 185 Bibliografie 192 Dicţionar de termeni folosiţi in lucrare 194 Bun de tipar 1 12 1971 Coli <le tipar 12,25 C Z 519 14 seria BAZELE MATEMATICE ALE CERCETĂRII OPERAŢIONALE